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The spontaneous breaking of the continuous symmetries of the two-dimensional (2D) electron-hole 
systems in a strong perpendicular magnetic field leads to the formation of new ground states and 
determines the energy spectra of the collective elementary excitations appearing over these ground 
states. In this review the main attention is given to the electron-hole systems forming coplanar 
magnetoexcitons in the Bose-Einstein condensation (BEC) ground state with the wave vector k = 0, 
taking into account the excited Landau levels, when the exciton-type elementary excitations coexist 
with the plasmon-type oscillations. At the same time properties of the two-dimensional electron 
gas (2DEG) spatially separated as in the case of double quantum wells (DQWs) from the 2D hole 
gas under conditions of the fractional quantum Hall effect (FQHE) are of great interest because 
they can influence the quantum states of the coplanar magnetoexcitons when the distance between 
the DQW layers diminishes. We also consider in this review the bilayer electron systems under 
conditions of the FQHE with the one half filling factor for each layer and with the total filling factor 
for two layers equal to unity because the coherence between the electron states in two layers is 
equivalent to the formation of the quantum Hall excitons (QHExs) in a coherent macroscopic state. 
This makes it possible to compare the energy spectrum of the collective elementary excitations of 
the Bose-Einstein condensed QHExs and coplanar magnetoexcitons. The breaking of the global 
gauge symmetry as well as of the continuous rotational symmetries leads to the formation of the 
gapless Nambu-Goldstone (NG) modes while the breaking of the local gauge symmetry gives rise 
to the Higgs phenomenon characterized by the gapped branches of the energy spectrum. These 
phenomena are equivalent to the emergence of massless and of massive particles, correspondingly, 
in the relativistic physics. The application of the Nielsen-Chadha theorem establishing the number 
of the NG modes depending of the number of the broken symmetry operators and the elucidation 
when the quasi-NG modes appear are demonstrated using as an example related with the BEC of 
spinor atoms in an optical trap. They have the final aim to better understand the results obtained in 
the case of the coplanar Bose-Einstein condensed magnetoexcitons. The Higgs phenomenon results 
in the emergence of the composite particles under the conditions of the FQHE. Their description 
in terms of the Ginzburg-Landau theory is remembered. The formation of the high density 2D 
magnetoexcitons and magnetoexciton-polaritons with point quantum vortices attached is suggested. 
The conditions in which the spontaneous coherence could appear in a system of indirect excitons 
in a double quantum well structures are discussed. The experimental attempts to achieve these 
conditions, the main results and the accumulated knowledge are reviewed. 

PACS numbers: 71.35.Lk, 67.85.Jk 
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I. INTRODUCTION 

The collective elementary excitations of the two- 
dimensional (2D) electron-hole (e-h) systems in a strong 
perpendicular magnetic field are discussed in the frame 
of the Bogoliubov theory of quasiaverages [1] taking into 
account the phenomena related with the spontaneous 
breaking of the continuous symmetries. The main re- 
sults in this field have been obtained thanks to the fun- 
damental papers by Goldstone [2], Nambu [3], Higgs [4] 
and Weinberg [5]. These investigations were influenced 
by the success of the theory of superconductivity devel- 
oped originally by Bardin, Cooper and Schriffer [6], re- 
fined later by Bogoliubov [1] as well as by the micro- 
scopic theory of superfluidity proposed by Bogoliubov 
[1]. The specific implementation of these concepts and 
theorems in the case of 2D magnetoexcitons with direct 
implication of the plasmon-type excitations side-by-side 
with the exciton-type branches of the energy spectrum is 
the main topic of the present review. The coplanar elec- 
trons and holes in a strong perpendicular magnetic field 
at low temperatures form the magnetoexcitons, when the 
Coulomb interaction between electrons and holes lying 
on the lowest Landau levels (LLLs) plays the main role. 
However, when the electrons and holes are spatially sep- 
arated on the different layers of the double quantum well 
(DQW) the Coulomb e-h interaction diminishes, and the 
two-dimensional electron gas (2DEG) on one layer and 
the two-dimensional hole gas (2DHG) on another layer 
are formed. Their properties under the conditions of the 
fractional quantum Hall effect (FQHE) can influence the 
properties of the 2D magnetoexcitons. To the best of our 
knowledge these aspects of the magnetoexciton physics 
were not discussed in literature. A short review is given 
on the Bose-Einstein Condensation (BEC) of the quan- 
tum Hall excitons (QHExs) arising in the bilayer electron 
systems under the conditions of the FQHE at one half fill- 
ing factor for each layer and the total filling factor equal 
to unity for both layers. This enables us to compare 
the phenomenon of the BEC of coplanar magnetoexci- 
tons and of QHExs. Such comparison provides better 
understanding of the underlying physics and allows to 
verify accuracy of the made approximations. Because the 
point vortices play an important role in the understand- 
ing of the FQHE the corresponding additional informa- 
tion should be included. The possibility to consider the 
BEC at T = as an estimate for the finite temperatures 
below the Berezinskii-Kosterlitz-Thouless phase transi- 
tion is suggested. The article is organized as follow. In 
Section 2 the Bogoliubov theory of the quasiaverages is 
overviewed. Section 3 is devoted to the Goldstone the- 
orem. The Nambu- Goldstone modes arising under the 
condition of BEC of the sodium atoms are enumerated 



in Section 4. The breaking of the local gauge symme- 
try and the Higgs phenomenon are discussed in Section 
5. Section 6 is devoted to the quasi-Nambu- Goldstone 
modes. In Section 7 the Ginzburg-Landau theory for the 
FQHE is formulated. The 2D point quantum vortices 
are described in Section 8. The existence of the statis- 
tical gauge vector potential generated by the vortices is 
considered in Section 9. The BEC of QHExs and the en- 
ergy spectrum of elementary excitations under these con- 
ditions are discussed in Section 10. Section II contains 
the main results concerning the energy spectrum of the 
exciton and plasmon branches of the collective elemen- 
tary excitations of the Bose-Einstein condensed coplanar 
magnetoexcitons. 



II. BOGOLIUBOV'S THEORY OF 
QUASIAVERAGES 

N.N. Bogoliubov [1] has demonstrated his concept of 
quasiaverages using the ideal Bose-gas model with the 
Hamiltonian 
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here at, are the Bose operators of creation and anni- 
hilation of particles, and \i is their chemical potential. 
The occupation numbers of the particles are 
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where fi < and /3 = 1/kT. 

In the normal state, the density of particles in 
the thermodynamic limit at ji = becomes n = 

2.612 {mk B Tf /2 /(2irh 2 f /2 . At this point, the Bose- 
Einstein condensation occurs and a finite value of the 
density of condensed particles appears in the thermody- 
namic limit 



no = lim — -; /i = — kgTln 
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The operators aj and ao asymptotically become 
c— numbers, when their commutator 
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asymptotically tends to zero and their product is equal 
to uq. Then one can write 
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On the other hand, the regular averages of the operators 
a' and a in the Hamiltonian (1) are exactly equal to 
zero. It is the consequence of the commutativity of the 
operator H and the operator of the total particle number 
N as follows 

# = 5>fc°k; [H,N]=0. (6) 

k 

As a result, the operators H is invariant with respect to 
the unitary transformation 

U = (7) 

with an arbitrary angle (p. This invariance is called 
gradient invariance of the first kind or gauge invari- 
ance. When 4> does not depend on the coordinate x, 
we have the global gauge invariance and in the case <f>(x) 
it is called local gauge invariance [2-8] or gauge invari- 
ance of the second kind. The invariance (7) implies 
H = U^HU; U^a U = e^a , which leads to the fol- 
lowing average value 

(a ) = Tr (a e- lm ) = Tr (a Ue- l)H U^) = 
= Tr(U^a Q Ue-P H ) = e* (a ) ; 
(1 - e*) (oo> = 

Because (p is an arbitrary angle, there are the selection 
rules: 

(a )=0; (4)=0 (8) 

The regular average (8) can also be obtained from the 
asymptotical expressions (5) if they are integrated over 
the angle a. This apparent contradiction can be resolved 
if Hamiltonian (1) is supplemented by additional term 

- v (a\e^ + aoe - * 1 ) W, v > 0, (9) 

where ip is the fixed angle and v is infinitesimal value. 
New Hamiltonian has the form 

(fc2 h.2 \ 
— - A a\a k - v (a\e^ + aoe"**) W. 

(10) 

It does not conserve the condensate number. Now the 
regular average values of the operators a\ and ao over 
the Hamiltonian H v ^ differ from zero, i.e., (ao) Hlp ^ 

and ^ a o^ 7^ 0. The definition of the quasiaverages 

designated by < a > is the limit of the regular average 
{ao) Hip when v tends to zero 

< ao >= lim (a ) Hu ip . (11) 

It is important to emphasize that the limit v — > must be 
effectuated after the thermodynamic limit V — > oo, A^o 
oo. In the thermodynamic limit, /i is also infinitesimal, 



and it is possible to choose the ratio of two infinitesimal 
values [i and v obtaining a finite value 

- - = y/r*. (12) 
M 

To calculate the regular average (ao) Hv one needs to 
represent the Hamiltonian (10) H v ^ in a diagonal form 
with the aid of the canonical transformation over the 
amplitudes 

a = --e^W + a Q] (13) 
A 4 

a k = a k ; k ^ 0. 

In terms of the new variables the Hamiltonian H v ^ takes 
the form 

t ^fh 2 k 2 \ t u 2 v „ A . 

H v>4> = -^o^ + 2^y-^--tija , i a % +— . (14) 

In the diagonal representation (14), the regular average 
value (ao)^^ exactly equals to zero, while the value 
(cio)h 1 , v i s equal to the first term on the right-hand side 
of formulas (13). 

As a result, the quasiaverage < a > is 

< a >= lim (ao) H - ^/W e^ (15) 

It depends on the fixed angle ip and does not depend on 
v. The spontaneous global gauge symmetry breaking is 
implied when the phase ip of the condensate amplitude 
in Hamiltonian (10) is fixed. 

When the interaction between the particles is taken 
into account, these differences appear for other ampli- 
tudes as well. They give rise to the renormalization of 
the energy spectrum of the collective elementary excita- 
tions. In such a way, the canonical transformation 

a k = 7Ao4,oe 4V + a k (16) 

introduced for the first time by Bogoliubov [1] in his the- 
ory of superfluidity, has a quantum-statistical foundation 
within the framework of the quasiaverage concept. At 
T = the quasiaverage < a > coincides with the aver- 
age over the quantum-mechanical ground state, which is 
the coherent macroscopic state [9] . 

The phenomena related to the spontaneous breaking 
of the continuous symmetry play an important role in 
statistical physics. 

Some elements of this concept, such as the coherent 
macroscopic state with a given fixed phase and the dis- 
placement canonical transformation of the field operator 
describing the Bose-Einstein condensate, were introduced 
by Bogoliubov in the microscopical theory of superfluid- 
ity [1] and were generalized in his theory of quasiaverages 
[1] noted above. 

The brief review of the gauge symmetries, their sponta- 
neous breaking, Goldstone and Higgs effects will be pre- 
sented below following the Ryder's monograph [7] and 
Berestetskii's lectures [8]. 
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III. GOLDSTONE'S THEOREM 

Goldstone has considered a simple model of the com- 
plex scalar Bose field to demonstrate his main idea. In 
the classical description the Lagrangian is 

-*<*•«• <■"> 

The potential energy V(fa has the form 

V(fa = m 2 fa<f) + \(fafa 2 ; A>0, (18) 

where m 2 is considered as a parameter only, rather than a 
mass term, A is the parameter of self-interaction, whereas 
the denotations and x M mean 

x M = (ct, x); x^ = (ct,—x); (19) 

The Lagrangian is invariant under the global gauge trans- 
formation 

<t> = e iK <t>; L(<j>) = L{4>'); A-const. (20) 

It has the global gauge symmetry. The ground state is 
obtained by minimizing the potential as follows 

^M =m V + 2A^| 2 . (21) 

Of interest is the case m 2 < 0, when the minima are 
situated along the ring 

H 2 = -^ = « 2 ; 101= a; a > 0. (22) 

The function V(<f>) is shown in Fig. 1 being plotted 
against two real components of the fields fa and fa. 
There is a set of degenerate vacuums related to each other 
by rotation. The complex scalar field can be expressed 
in terms of two scalar real fields, such as p(x) and 9(x), 
in polar coordinates representation or in the Cartesian 
decomposition as follows 

fax) = p(x)e ie ^ = (fa( x ) + ifa( x )) _L (23) 

The Bogoliubov-type canonical transformation breaking 
the global gauge symmetry is 

M=a + * {x) ^ X) =Wx) + a)e«V (24) 

The new particular vacuum state has the average (<f>) = 
a with the particular vanishing vacuum expectation val- 
ues (4>i)q = (4>2)o = (p')o = (^')o = ^- ^ means the 
selection of one vacuum state with infinitesimal phase 
9' — > 0. As was pointed in [7], the physical fields are the 
excitations above the vacuum. They can be realized by 
performing perturbations about \<f>\ — a. Expanding the 



V 

t 




FIG. 1. The potential V(<f>) with the minima at \<j>\ — a and 
a local maximum at <f> — 

Lagrangian (17) in series of the infinitesimal perturba- 
tions 9', p', (j)' 2 and neglecting by the constant terms, 
we obtain 

L = \ (Vi) (d^[) + \ (d^' 2 ) - 2\a 2 <j> 12 - 

-V2A0' 1 (0' I + 2 2 )-^(0' 1 2 + 2 2 ) 2 (25) 

or in polar description 

L = ( V) (d»p f ) + (p 1 + af (d,9') (d»6') 

- [\p lA + 4aAp' 3 + 4\a 2 p 12 - Aa 4 ] (26) 

Neglecting by the cubic and quartic terms, we will see 
that there are the quadratic terms only of the type 
4Xa 2 p' 2 and 2Xa 2 <p' 1 2 , but there is no quadratic terms 
proportional to 9' 2 and (j>' 2 2 ■ For real physical problems, 
for example, for the field theory, the field components 
f/4 and p' represent massive particles and dispersion laws 
with energy gap, whereas the field components tf>' 2 and 9' 
represent the massless particles and gapless energy spec- 
trum. 

The main Goldstone results can be formulated as follows 
m 2 , = 4Aa 2 , m|/ = 2Aa 2 ; 

m 2 , -0; m 2 , = 0; (27) 

The spontaneous breaking of the global gauge symmetry 
takes place due to the influence of the quantum fluctua- 
tions. They transform the initial field <j> with two mas- 
sive real components 4>i and fa, and a degenerate ground 
state with the minima forming a ring into another field 
with one massive and other massless components, the 
ground state of which has a well defined phase without 
initial symmetry. 

The elementary excitations above the new ground state 
changing the value (p) = a are massive. It costs energy 
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to displace p 1 against the restoring forces of the potential 
V(p). But there are no restoring forces corresponding to 
displacements along the circular valley \(f>\ = a formed by 
initial degenerate vacuums. 

Hence, for angular excitations 6'oi wavelength, A we have 
to <~ A -1 — > as A — > oo. The dispersion law is uo <~ ck 
and the particles are massless [7]. The 8' particles are 
known as the Goldstone bosons. This phenomenon is 
general and takes place in any order of perturbation the- 
ory. The spontaneous breaking of a continuous symmetry 
not only of the type as a global gauge symmetry but also 
of the type of rotational symmetry entails the existence 
of massless particles referred to as Goldstone particles 
or Nambu-Golstone gaplcss modes. This statement is 
known as Goldstone theorem. It establishes that there 
exists a gapless excitation mode when a continuous sym- 
metry is spontaneously broken. The angular excitations 
9' are analogous to the spin waves. The latter represent 
a slow spatial variation of the direction of magnetization 
without changing of its absolute value. Since the forces in 
a ferromagnetic are of short range, it requires a very lit- 
tle energy to excite this ground state. So, the frequency 
of the spin waves has the dispersion law to — ck. As 
was mentioned by Ryder [7], this argument breaks down 
if there are long-range forces like, for example, the 1/r 
Coulomb force. In this case, we deal with the maxwellian 
gauge field with local depending on x gauge symmetry 
instead of global gauge symmetry considered above. 

After the specific application of the above statement 
will be demonstrated following Refs. [10-16], where 
the spinor Bose-Einstein condensates were discussed, we 
will consider the case of Goldstone field <j) and of the 
maxwellian field with local gauge symmetry. 



IV. BOGOLIUBOV'S EXCITATIONS AND THE 
NAMBU-GOLDSTONE MODES 

The above formulated theorems can be illustrated us- 
ing the specific example of the Bose-Einstein condensed 
sodium atoms 23 Na in an optical-dipole trap following 
the investigations of Murata, Saito and Ueda [10] on 
the one side and of Uchino, Kobayashi and Ueda [11] 
on the other side. There are numerous publications on 
this subject among which should be mentioned [12-17]. 
The sodium atom 23 Na has spin / = 1 of the hypcrfmc 
interaction and obey the Bose statistics. Resultant spin 
of the interacting bosons with / = 1 is F which takes 
the values F — 0,1,2. The contact hard-core interac- 
tion constant gF — ^Trh 2 ap/M is characterized by s-wave 
scattering length a^, which is not zero for F ~ when 
two atomic spins form a singlet, and for F — 2, when 
they form a quintuplet. The constant go and gi enter 
the combinations cq = (go + 1g-i)j2> and c\ = (<?2 — <?o)/3 
which determine the Hamiltonian. The description of 
the atomic Bose gas in an optical-dipole trap is possible 
in the plane- wave representation due to the homogeneity 
and the translational symmetry of the system. It means 



that the components of the Bose field operator ip m (r) can 
be represented in the form: 

V> m (r) = -^^a fcm e 4fc> (28) 

where V is the volume of the system, and a^ m is the 
annihilation operator with the wave vector k and the 
magnetic quantum number m, which in the case / = 1 
takes three values 1, 0, -1. The spinor Bose-Einstein con- 
densates were realized experimentally by the MIT group 
[12] for different spin combinations using the sodium 
atoms 23 Na in a hyperfine spin states \ f — l,m,f = — 1) 
in a magnetic trap and then transforming them to the 
optical-dipole trap formed by the single infrared laser. 
The Bose-Einstein condensates were found to be long- 
lived. Some arguments concerning the metastable long- 
lived states were formulated. The states may appear if 
the energy barriers exist, which prevent the system from 
direct evolving toward its ground states. If the thermal 
energy needed to overcome these barriers is not available, 
the metastable state may be long-lived and these events 
are commonly encountered. Even the Bose-Einstein con- 
densates in the dilute atomic gases can also be formed 
due to the metastability. Moreover, in the gases with at- 
tractive interactions the Bose-Einstein condensates may 
be metastable against the collapse just due to the en- 
ergy barriers [12]. Bellow we will discuss the Bogoliubov- 
type collective elementary excitations arising over the 
metastable long-lived ground states of the spinor-type 
Bose-Einstein condensates (BEC-tes) following Ref.[10, 
11], so as to demonstrate the formation of the Nambu- 
Goldstonc modes. 

The Hamiltonian considered in [10] is given by formu- 
las (3) and (4), and has the form 

H = E ( £ k - P m + l™ 2 ) a l m a km + 
k,m 

^E^+^E^ (^) 

k k 

Here the following designations were used 
H 2 k 2 

Ek = I'M - ' C ° = ^° + 2 - 92 )/ 3 ' Cl = (- 92 ~ 5 °)/ 3 

h = E 4,^: ?= If*,?, /*); (30) 

q,m 

q,'m,n 
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The repeated indexes mean summation over 1,0,-1. The 
symbol :: denotes the normal ordering of the operators. 
The coefficient p is the sum of the linear Zeeman energy 
and of the Lagrangian multiplier, which is introduced to 
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set the total magnetization in the z direction to a pre- 
scribed value. This magnetization is conserved due to 
the axisymmetry of the system in a magnetic field, q 
is the quadratic Zeeman effect energy, which is positive 
in the case of spin / = 1 for 23 Na and 87 Rb atoms. The 
spin-spin interaction is of ferromagnetic-type with c\ < 
for / = 1 87 Rb atoms and is antiferromagnetic-type with 
ci > for / = 1 23 Na atoms [10]. Taking into account 
that in many experimental situations the linear Zeeman 
effect can be ignored and the quadratic Zeeman effect 
term q can be manipulated experimentally, in [11] both 
cases of positive and negative q at p = were investi- 
gated for spin —1 and spin —2 Bosc-Einstein condensates 
(BECs). We restrict ourselves to review some spinor 
phases with spin —1 discussed in [11] so as to demonstrate 
the relations between the Nambu-Goldstone(NG) modes 
of the Bogoliubov energy spectra and the spontaneous 
breaking of the continuous symmetries. The description 
of the excitations is presented in Refs. [10, 11] in the 
number-conserving variant of the Bogoliubov theory [1]. 
There is no need to introduce the chemical potential as 
a Lagrangian multiplier in order to adjust the particle 
number to a prescribed value. The BEC takes place on 
a superposition state involving the single-particle states 
with wave vector k = and different magnetic quantum 
numbers 

lo = £&»°Lh' ac >; £iu 2 = i (si) 

m m 

The order parameter has a vector form and consists of 
three components: £ = (£i,£o>£-i)- The vacuum state 
vac) means the absence of the atoms. The ground state 
wave function of the BEC-ed atoms is given by the for- 
mula (8) of Ref. [11] 

^ a) = ~m I ^ &»°o,m] \ vac ) ( 32 ) 



In the mean- field approximation the operators a m , a m 
are replaced by the c— numbers £ m \/No, where No is the 
number of the condensed atoms. After this substitution, 
the initial Hamiltonian loses its global gauge symmetry 
and does not commute with the operator N. The order 
parameters £ m are chosen so that to minimize the expec- 
tation value of the new Hamiltonian and its ground state 
and to satisfy the normalization condition ^ | £ m | 2 = 1. 

m 

To keep the order parameter of each phase unchanged 
it is necessary to specify the combination of the gauge 
transformation and spin rotations [11]. This program 
was carried out in [18-21]. 

The initial Hamiltonian (29) in the absence of the ex- 
ternal magnetic field has the symmetry U(l) x SO (3) 
representing the global gauge symmetry U(l) and the 
spin- rotation symmetry 5*0(3). The generators of these 
symmetries are referred to as symmetry generators and 



have the form 

N = J dxft m {x)j) m { X ) =J2 a i, m a k,m 

k,m 

= J dxi> m (x)fiJ n (x); j =x,y,z (33) 

Unlike the 5*0(3) symmetry group with three generators 
F x , F v and F z , the 50(2) symmetry group has only one 
generator F z which describes the spin rotation around 
the z axis and looks as follows: 

F z = mat at (34) 

k.m 

In the presence of an external magnetic field, the sym- 
metry of the Hamiltonian is U(l) x 50(2). The break- 
ing of the continuous symmetry means the breaking of 
their generators. The number of the broken generators 
(BG) is denoted as Nbg- There are 4 generators in the 
case of U(l) x 50(3) symmetry and two in the case of 
U(l) x 50(2)symmetry. 

The phase transition of the spinor Bose gas from the 
normal state to the Bosc-Einstein condensed state was 
introduced mathematically into Hamiltonian (29) using 
the Bogoliubov displacement canonical transformation, 
when the single-particle creation and annihilation oper- 
ators with the given wave vector k, for example, k = 0, 
were substituted by the macroscopically c— numbers de- 
scribing the condensate formation. The different super- 
positions of the single-particle states determine the struc- 
ture of the finally established spinor phases [11]. Nielsen 
and Chadha [17] formulated the theorem which estab- 
lishes the relation between the number of the Nambu- 
Goldstonc modes, which must be present between the 
amount of the collective elementary excitations, which 
appear over the ground state of the system if it is formed 
as a result of the spontaneous breaking of the Nbg con- 
tinuous symmetries. The number of NG modes of the 
first type with linear (odd) dispersion law in the limit of 
long wavelengths denoted as Nj being accounted once, 
and the number Nji of the NG modes of the second type 
with quadratic (even) dispersion law at small wave vec- 
tors, being accounted twice give rise to the expression 
Ni + 2Nn, which is equal to or greater than the number 
Nbg °f the broken symmetry generators. The theorem 
[17] says 

Ni + 2N U > N BG (35) 

The theorem has been verified in [11] for multiple ex- 
amples of the spin —1 and spin —2 Bose-Einstein con- 
densate phases. In the case of spin —2 nematic phases, 
the special Bogoliubov modes that have linear dispersion 
relation but do not belong to the NG modes were re- 
vealed. The Bogoliubov theory of the spin —1 and spin 
—2 Bose-Einstein condensates (BECs) in the presence of 
the quadratic Zeeman effect was developed by Uchino, 
Kobayashi and Ueda [11] taking into account the Lee, 
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Huang, Yang (LHY) corrections to the ground state en- 
ergy, pressure, sound velocity and quantum depletion of 
the condensate. Many phases that can be realized ex- 
perimentally were discussed to examine their stability 
against the quantum fluctuations and the quadratic Zee- 
man effect. The relations between the numbers of the 
NG modes and of the broken symmetry generators were 
verified. A brief review of the results concerning the spin 
— 1 phases of [11] is presented below so as to demonstrate, 
using these examples, the relations between the Bogoli- 
ubov excitations and the Nambu-Goldstone modes. 

The first example is the ferromagnetic phase with c\ < 
0, q < and the vector order parameter 

f = (1,0,0) (36) 

The modes with m — and m = — 1 are already diag- 
onalizcd, whereas the mode m = 1 is diagonalized by 
the standard Bogoliubov transformation. The Bogoli- 
ubov spectrum is given by formulas (33) and (34) of Ref. 
[11] 

E k,i = ^kiSk + Mco + C!)); (37) 

E k,o = £ k-^ E k,-i = £ k- 2c ^ n 

The Ej: 1 mode is masslcss. In the absence of a magnetic 
field, when q = 0, the mode m = is also massless with 
the quadratic dispersion law. The initial symmetry of the 
Hamiltonian before the phase transition is U(l) x 50(3), 
whereas the final, remaining symmetry after the process 
of BEC is the symmetry of the ferromagnetic i.e. 50(2). 
From the four initial symmetry generators N, F x , F v and 
F z remains only the generator F z of the 50(2) symme- 
try. The generators F x and F y were broken by the ferro- 
magnet phase, whereas the gauge symmetry operator N 
was broken by the Bogoliubov displacement transforma- 
tion. The number of the broken generators N, F x , F v is 
three, i.e., Nbg = 3. In this case Ni = 1, Nu — 1 and 
Nj + 2Njj = 3, being equal to Nbg = 3. The equality 
Nj + 2Nji = Nbg takes place. In the presence of an 
external magnetic field, with q ^ 0, the initial symme- 
try before the phase transition is U(\) x 50(2) with two 
generators N and F z , whereas after the BEC and the 
ferromagnetic phase formation the remained symmetry 
is 50(2). Only one symmetry generator N was broken. 
It means Nbg = 1, iVj = 1 and Njj = 0. The equality 
Nj + 2Njj = Nbg also takes place. 

The condition (co + a) > to be hold is required for 
m = 1 the Bogoliubov mode to be stable. It ensures the 
mechanical stability of the mean-field ground state. Oth- 
erwise, the compressibility would not be positive definite 
and the system would become unstable against collapse. 
In the case q > 0, c\ > and (co + c\) < the state 
would undergo the Landau instability for the m = and 
m = — 1 modes with quadratic spectra and the dynami- 
cal instability for the m = 1 mode with a linear spectrum 
(36) of Ref. [11]. 



There are two polar phases. One with the parameters 

I* 3 = (0,1,0); q > 0; q + 2nci > (38) 
and the other with the parameters 

r' = -J= (1,0,1); q<0; Ci > (39) 

These two polar phases have two spinor configurations 
which are degenerate at q = and connect other by 
U(l) x 50(3) transformation. However, for nonzero q 
the degeneracy is lifted and they should be considered 
as different phases. This is because the phase P has 
a remaining symmetry 50(2), whereas the phase P' is 
not invariant under any continuous transformation. The 
number of NG modes is different in each phase and the 
low-energy behavior is also different. Following formu- 
las (40)- (42) of [11] the density fluctuation operator akd 
and the spin fluctuation operators dkj x and dkj y were 
introduced 

1 / > 

akd — flfe,o; 0-kJ x — ~^( ak ' 1 + a <^-l); 

a-hjy = y|( a M _ °fc -i); ( 40 ) 
Their Bogoliubov energy spectra are 

E u = V^( £ * + 2c °); ( 41 ) 

E kj 3 = V /(£fc + g)(^ + 9 + 2n Cl ) ; 

In the presence of an external magnetic field, the initial 
symmetry is U(l) x 50(2), whereas after the BEC and 
the formation of the phase P with g ^ the remaining 
symmetry is also SO (2). Only the symmetry U(l ) and 
its generator N were broken during the phase transition. 
It means we have in this case Nbg — 1) Nj — 1 and 
N H = 0. The equality Ni + 2N H = N B g holds. Den- 
sity mode is massless because the U(l) gauge symmetry 
is spontaneously broken in the mean- field ground state, 
while the transverse magnetization modes f x and f v are 
massive for non zero q, since the rotational degenera- 
cies about the x and y axes do not exist being lifted by 
the external magnetic field. In the limit of infinitesimal 
q — > nevertheless nonzero, the transverse magnetization 
modes f x and f y become massless. It occurs because be- 
fore the BEC in the absence of an external magnetic field 
the symmetry of the spinor Bose gas is U(l) x 50(3), 
whereas after the phase transition it can be considered 
as a remaining symmetry 50(2). The generators N, F x , 
E v were broken, whereas the generator F z remained. In 
this case we have Nbg = 3, Nj = 3 and Njj = the 
equality looks as 3 = 3. 

In the polar phase P 1 with the parameters (39) the 
density and spin fluctuation operators were introduced 
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by formulas (57)-(59) of Rcf. [11] 

1 , , 

flfed = -^={a,k,i + a k,-i); a,kj x = »fc,o; 

a kj y = -j={a k s - a k -i); 



(42) 



with the Bogoliubov energy spectra described by formu- 
las (65)-(67) [11]: 



E k,d = \/ £ k( £ k + 2nC ^ E kJ, = \/ £ k( £ k + 2nC ^ 



E kJ x = \ ( £ k~l)( £ k~1 + 2nc ^ 



(43) 



At q < in contrast to the case q > one of the spin 
fluctuation mode E^ j. becomes massless. The initial 
symmetry of the system is U (1) x 50(2). It has the sym- 
metry generators N and F z . They are completely broken 
during the phase transition. After the phase transition 
and the P' phase formation there are not any symmetry 
generators. The number of the broken generator is 2 ( 
Nbg = 2), whereas the numbers Nj and Nu are 2 and 
0, respectively. As in the previous cases, the equality 
occurs in the Nielsen and Chadha rule. For the Bogoli- 
ubov spectra to be real the condition q < 0, c > and 
c\ > must be satisfied, otherwise, the state £ p will be 
dynamically unstable. 

Side by side with the spinor-type three-dimensional 
(3D) atomic Bosc-Einstein condensates in the optical 
traps, we will discuss also the case of the Bose-Einstein 
condensation of the two-dimensional (2D) magnetoexci- 
tons in semiconductors [22-25]. The collective elemen- 
tary excitations under these conditions were investigated 
in [26-31] and will be described in Section 11. As was 
shown above, the spontaneous symmetry breaking yields 
Nambu-Goldstone modes, which play a crucial role in 
determining low-energy behavior of various systems [5, 
32-38] . Side by side with the global gauge symmetry the 
local symmetry does exist. 



V. SPONTANEOUS BREAKING OF THE 
LOCAL GAUGE SYMMETRY AND THE HIGGS 
PHENOMENON 

The interaction of the electrons with the electromag- 
netic field can be described introducing into the La- 
grangian the kinetic momentum operators instead of 
canonical ones what is equivalent to introduce the co- 
variant derivatives D instead of the differential ones d. 
They are determined in Ref. [8] as 

,1 d 



D = d--A;A=(<p,A) 



(44) 



where ip and A are the scalar and vector potentials of the 
electromagnetic field (EMF). Below we will use also the 



denotations of Ref. [7] 

x" = (ct,r); x M = (ci,-r); A^ = (<p,-A); A" = {<p,A); 



a - JL 



1 d 



9/19 



cmV'^-wAiw-^ (45) 



e : 

c 



E 

= I —,-P 



i_ a 2 

c 2 dt' \ c ~ ) ~ ' \ c 

The Lagrangian of the free EMF has the form [7] 

1 

4 



Lemf — ——F IIV F ,1V 



(46) 



being expressed through the antisymmetric tensors F^ 
and F^ . They are determined as four-dimensional curls 



of A M and A 11 . 



1 ILV 



dpAv-dvAn; F^ = d»A v -d v A» (47) 



The full Lagrangian of the electrons and EMF reads [7] 

ie 



- m 



2i* 



—A 



he 



.A" 



- (48) 



As before m 2 is a parameter so that in the case m 2 < 
and in the absence of the EMF vacuum values are deter- 
mined by the formula (22). 

The invariance of the Lagrangian (48) under the trans- 
formation [8] 



(f>'(x) = (f)(x)e 



i6(x) 



(49) 



in the presence of the EMF can be achieved only under 
the concomitant transformation of its potential in the 
form [8] 



he , 



A\x) = A(x) + —d9(x) 



(50) 



Indeed in this case: the Lagrangian (48) remains invariant 
[7] as follows 



F' -- 



ie 
he 
ie 
he 



A' 1 
A'" 



(51) 



Introducing the gauge transformation of the field func- 
tion (24) and expanding the Lagrangian in power series 
on the small physical fields <$>\ and 4>' 2 we obtain the con- 
stant, quadratic, cubic and quartic terms. The quadratic 
part looks as [7] 



+ 1 -(dM 2 + 1 -(d^' 2 f- 

-2\a 2 <i>' 2 + V2eaA^d^' 2 



(52) 
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The second term is proportional to A^A^. It indicates 
that the photon becomes massive. The scalar field <ft x is 
also a massive one. The field <j)' 2 takes part in the mixed 
term A^d^' 2 and can be eliminated by the supplemen- 
tary gauge transformation (50). Following Ref [7] the 
Lagrangian (52) can be presented in the form 



L 2 = -- A F^F^ + e 2 a 2 A^ + 
+ \ M) 2 - 2\a 2 cj>' 2 



(53) 



It contains two fields only: the photon with longitudinal 
component and spin 1 and field 4> x with spin 0. They are 
both massive. The field </>' 2 , which in the case of sponta- 
neous breaking of the global symmetry became massless 
forming a Goldstone boson, in this case disappeared. The 
photon became massive. This phenomenon is called the 
Higgs phenomenon [7]. 

One possible illustration of the described above effect 
will be considered below following the paper by Halperin, 
Lee and Read [39] . They considered the two-dimensional 
(2D) system of spinless electrons under the conditions of 
the quantum Hall effect. Then the Hamiltonian H = 
K + V consists of the kinetic energy operator K 



K = 



1 



2m e 



<P~?ft e C?)[-iK$ + -l(^)] 2 ^ e (^) (54) 



with 2D electrons with the mass m e and the charge — e 
situated in a uniform external perpendicular magnetic 
field B with the vector potential A(r). The potential 
energy operator V depends on the Coulomb interaction 
between the electrons. The creation and annihilation op- 
erators iplif), ■0 e (r)obey to the Fermi statistics as was the 
case of Ref. [39] , but we will consider following the Ref. [40] 
a more general case including also the Bose statistics 

± = <5 2 (^ - f ) (55) 

The signs ± correspond to the Fermi and Bose statistics. 
In Ref. [39] the new " quasiparticle" operators f/^(r), tp(r) 
were introduced by the relations 

= ^(^)e" imS ;V^) = e lmQ ^(^) (56) 
with an integer number m and with the phase operator 

u)(r) = J d 2 r l 9(r- f f )p{f f ) (57) 

It depends on the angle 9(r — f) between the vector f — f" 
and the in-plane axis x being determined by the formula 



0(f— f) = arctan 



y-y 



(58) 



x — x' 

and by the density operator p{f) 

i>(f) = ^(r')i(fO = ft(?)4>(?) (59) 



These operators have the properties 

[TP e (r),p(r')]=TP e (r)5 2 (!>-7) 

$ e (f),w(?)] = &(»W-0 (60) 
4, e (r)LJ n (7) = (w(^) + 9(f- 7)) n 4> e (7) 
V> (■r)e imC ( ? '- ) = e "" e (^ ? + ? ') e im2 ( ? ')^, (f) 

It will be shown below that m is the number of point 
vortices attached to each bare initial particle forming to- 
gether with it a composite particle (CP). The statistics of 
the CPs depends on the statistics of the initial particles 
and on the number m of the attached vortices. Finally, 
we will calculate the commutators of the operators if)' (r) , 
i>(r) with the requirement that it will be S 2 (f— r') as fol- 
lows: 

$(?),^(i*)]± 

(61) 



xe 



-im{uj{r )~ uJ(f)) 



= 8 2 {-? -7) 



Here we have taken into account the relation Qif 1 — r) — 
9{r — f") = 7r for f ^ f 1 . One can observe that the CPs 
represented by the operators ^(r), i>{r) are composite 
fcrmions (CFs) if the bare initial particles are fermions 
and the number of vortices m is even as well as in the case 
when the initial particles are bosons and the number of 
vortices m is odd. In the same way the CPs are composite 
bosons (CBs) if the initial particles are fermions and the 
number of vortices is odd, or if the initial particles are 
bosons and the number of vortices m is even. 
The kinetic energy operator K in terms of the operators 
tp^r) and ip(r) is 

h 2 f p - 

K= / d 2 r^(r)e tmC, ^[-iV+— 4(r)] 2 e-~' r ^(r) 

2m e J he 

(62) 

It can be transformed taking into account that 



(-iV + — A{r))e- lm ^ip(r) 

he 



= e - im ^ ( r^_ lV + j-A(?) - mVw(r))V(r) (63) 

what leads to the formula 

h 2 f e 
K = - — / d 2 f^(f)[-iV + —A(f) - raVw(r)]V(r) 

(64) 

It contains a supplementary vector potential a(r) named 
as statistical Chern-Simons gauge potential [41] deter- 
mined as 

d(f) = _!^vu;(f) = -— I d 2 ?V r 9(r- r*)p{?) 

(65) 
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Its calculation needs a special precaution as was pointed 



by Jackiw and Pi [41] because 



') is a multivalued 



function. They cautioned against the moving of V with 
respect to f out of the integral J d 2 r'd{r — f f )p(r f ), be- 
cause in general it is not correct. The integration cannot 
be interchanged with the differentiation. The reason for 
this is that the function 0(f '— f) is multivalued and the 
integration of 6(f — f) over the two-dimensional f plane 
requires specifying the cut in the space f , which begins 
at the point f. The range of the f integration depends 
on f and moving the f derivative outside of the f in- 
tegral gives an additional contribution. To avoid these 
complications the derivative V6(f— f ) is introduced into 
the integrand for the very beginning in the form 



V6>(f — f ) = —curl In |f ■ 



(66) 



As it was shown in Ref. [41] the curl of a scalar S in 
the 2D space is a vector and the curl of the vector a is a 
scalar as follows 



(curlSy = e ij djS; curia =E ij d^; i,j = 1,2 



12 



e 21 ;e n =e 22 =0 



(67) 



The statistical gauge vector potential a(r) can be tran- 
scibcd 



a(r) 



mhc 



J d 2 f'curlln 



?-f f \p{r f ) (68) 



what leads to the statistical gauge magnetic field b{r) 
b{r) = curla{r) = ^ J d 2 f" G ij G jk djd k In |f - ?\p(7) = 



mhc 



I 



d 2 r'A r ln\f-f'\p(r') (69) 



Taking into account the equality 

Alnr = 27T(5 2 (r) 

we obtain 



6(f) 



2nmhc 



(70) 



(71) 



Substituting the density operator /5(f) by its mean value 
with the fractional integer filling factor v equal 



2tt/ 2 



to v = 1/m with m > 1, and taking into account the 
he 

magnetic length I 2 — — determined by the external 

magnetic field B we will find the average value b and 
equality 



B + b = 



(72) 



what means that the resulting magnetic field is exactly 
zero. In this approximation the set of CPs does exist 
in zero magnetic field. If they are fermions their ground 



state will be a filled Fermi sea with the Fermi wave vector 
determined by the magnetic length. If they are bosons, 
they will undergo the BEC. 

In the Section 7 we will discuss the collective elemen- 
tary excitations above the ground state in the case of CBs 
on the base of Ginzburg-Landau theory. In Ref. [42] it was 
shown that applying the mean-field theory one must inte- 
grate out the short-distance fluctuations of the ip(r) field 
to obtain an effective action which describe the physics at 
distance scales larger than the magnetic length. It is sup- 
posed that the effective action is of the same form as the 
microscopic action, but with renormalized stiffness con- 
stant, bare mass and the effective interaction strength. 



VI. QUASI-NAMBU-GOLDSTONE MODES IN 
THE BOSE- EINSTEIN CONDENSATES 

The Goldstone theorem guarantees that the NG modes 
do not acquire mass at any order of quantum corrections. 
Nevertheless, sometimes soft modes appear, which are 
massless in the zeroth order but become massive due to 
quantum corrections. They were introduced by Weinberg 
[5], who showed that these modes emerge if the symme- 
try of an effective potential of the zeroth order is higher 
than that of the gauge symmetry and the idea was in- 
voked to account for the emergence of low-mass particles 
in relativistic physics. Following [32] now these modes are 
referred to as quasi-Nambu-Goldstone modes, in spite of 
the fact that their initial name introduced by Weinberg 
was pseudo-modes instead of quasi-modes. Georgi and 
Pais [33] demonstrated that the quasi-NG modes also oc- 
cur in cases in which the symmetry of the ground state 
is higher than that of the Hamiltonian [32] . This type of 
the quasi-Nambu-Goldstone modes is believed to appear, 
for example, in the weak-coupled limit of A phase of 3 He 
[37, 38]. 

The authors of [32] underlined that the spinor BEC 
are ideal systems to study the physics of the quasi-NG 
modes, because these systems have a great experimental 
manipulability and well established microscopic Hamil- 
tonian. It was shown in [32] that the quasi-NG modes 
appear in a spin-2 nematic phase. In the nematic conden- 
sate, three phases, each of which has a different symme- 
try, are energetically degenerate to the zeroth order [36] 
and the zeroth order solution has a rotational symmetry 
SO(5), whereas the Hamiltonian of the spin-2 condensate 
has a rotational symmetry SO(3). By applying the Bo- 
goliubov theory of the BEC under the assumption that 
the k = components of the field operators are macro- 
scopically occupied, it was shown that the order param- 
eter of the nematic phase has an additional parameter 
independent on the rotational symmetry. 

The ground state symmetry of the nematic phase at 
the zeroth order approximation is broken by quantum 
corrections, thereby making the quasi-NG modes mas- 
sive. The breaking of the SO (5) symmetry occurs. The 
number n of the quasi-NG modes was determined by 
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Georgi and Pais [33] in the form of the theorem. It was 
explained and represented in [32] as follows: 



order parameter <j> has the form 



n = dim(M) - dim(M) 



(73) 



where M is the surface on which the effective potential as- 
sumes its minimal values to the zeroth order and dim(M) 
is the dimension of this surface. The dimension dim(M) 
determines the number of the NG modes. This implies 
that M is a submanifold of M and n is the dimension of 
the complementary space of M inside M [32] . 

In the case considered by Goldstone, the dimension of 
the ring is 1 and the number of the NG modes is 1. This 
leads to the absence of the quasi-NG modes (n = 0). Re- 
turning to the case of 2D magnetoexcitons in the BEC 
state with small but nonzero wave vector k (k =/= 0) de- 
scribed by Hamiltonian (16) of [30], one should remember 
that both continuous symmetries existing in the initial 
form (10) [30] were lost. It happened due the presence of 
the term ij(dl + dg) in the frame of the Bogoliubov theory 
of quasiaverages. Nevertheless, the energy of the ground 
state as well as the self-energy parts Sjj(P, w), which 
determine the energy spectrum of the collective elemen- 
tary excitations depend only on the modulus of the wave 
vector k and do not depend at all on its direction. All 
these expressions have a rotational symmetry SO(2) in 
spite of the fact that Hamiltonian (16) of [30] has lost it. 
We believe that the condition described by Georgi and 
Pais [33] favoring the emergence of the quasi-NG modes. 
We are explaining the existence of the gapped, massive 
exciton-type branches of the collective elementary exci- 
tations obtained in our calculations just by these con- 
siderations. These questions will be discussed in Section 
11. 



VII. GINZBURG-LANDAU THEORY FOR THE 
FRACTIONAL QUANTUM HALL EFFECT 

In this section we will follow the collective mono- 
graph [43] dedicated to the fractional quantum Hall ef- 
fect (FQHE) , the clear and transparent candidatus scien- 
tiarum thesis by Enger [44] and many other papers cited 
below. The Landau theory of the second order phase 
transition [45] is based on the introduction of the order 
parameter, <p(rt) assuming that the free energy is a reg- 
ular function of <fr at least near the critical point. In the 
case of superconductors and superfluids the role of the 
order parameter is played by the condensate wave func- 
tions. The theory of superconductors was elaborated by 
Ginzburg and Landau [46] whereas for liquid helium by 
Ginzburg and Pitaevskii [47] . The microscopical founda- 
tions in the latter case were proposed by Pitaevskii [48] 
and by Gross [49] and can be found in the monograph 
by Nozieres and Pines [50] . The microscopical theory of 
superfluidity was firstly proposed by Bogoliubov as the 
model of weakly interacting Bose gas [51]. The density 
of the Hclmholz free energy /(r) expanded on the small 



f{r) = f + a\<f>\' 



H 4 



2m 



(74) 



In the case of superconductor it is necessary to include 
the effect of the applied electromagnetic field which can 
be done by substituting the canonical momentum p = 
—ihV by the kinetic momentum 



p- -A(r) 

c 



(75) 



where A is the vector potential, q is the charge of the 
Cooper pair, q = — 2e. 

The density of the Gibbs free energy including also the 
density of the magnetic field energy looks as 



5 (r)=./b + aH 2 + ^H 4 + 



1 

2m 



-ihv + —A 



r>2 



where B = rot A. Minimizing the total Gibbs energy 
G = J g(r)dr with respect to <j> and A gives 

J-(-ihv + — A) + c^ + /3|0| 2 = o 

2m \ cl 



'' -> ieti 4p 2 

-VxB = - (0*V0 - <t>V<f>*) - ^ 
Ho m rac z 



\<t>\ A (77) 



This is the Ginzburg-Landau equations, where /x is the 
magnetic permeability. 

The Ginzburg-Pitaevskii-Gross equation for the Bose- 
Einstein condensate wave function 0(r, t) is 



ih 



d <f>(r,t) 
dt 



-|^-?0(r,t)+A|0(r,i)| 2 0(r,t) (78) 



Separating the space and time parts <fi(r,t) — e~ lfJ,t (j)(r), 
and choosing the chemical potential fi = Apo, one can 
transform (78) into the equation 



(79) 



which is known as Gross-Pitaevskii equation or non- 
linear Schrodinger equation. As was mentioned in [51, 
52] the GL theory is needed also for the FQHE to better 
understand this phenomenon. 

The FQHE also is a remarkable example of the quan- 
tum effects observable on a macroscopic level similarly as 
superconductivity and superfluidity. All these phenom- 
ena have a ground state with non-zero density of particles 
and in all three cases there are quasiparticle excitations 
in the form of vortices. But there are some aspects of 
the FQHE, which are not present in the GL theories of 
superconductors and superfluids. First of them there is a 
gap in the spectrum of the collective elementary excita- 
tions, which leads to the incompressibility of the FQHE 
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systems. The second important difference is related with 
the properties of the vortices in the FQHE case. They 
play the role of the single-particle excitations and have 
finite creation energy, as opposed to the vortices in the 
superfluid He-II with an extensive creation energy of the 
vortex proportional to ln(i?/a), where R is the radius of 
the system and a is the vortex core. 
In addition the FQHE vortices have fractional charges 
[52] . In numerous papers some variants of the G-L the- 
ory for the FQHE were proposed starting with the La- 
grangian of the system containing the supplementary 
term known as Chern-Simons term. It describes the 
gauge vector potential generated by the vortices; which 
in their turn are induced by the flux quanta created by 
the external magnetic field B. Instead of Gibbs free en- 
ergy the action of the system is studied. 
Girvin [52], and Girvin and MacDonald [53] for the first 
time proposed a phenomenological variant of the GL the- 
ory writing the action S in the form 



d 2 r{ (-ihS7 + -A 1 (r)ip(r) 



where 



+i(V>*(r)V(r) - n Q )4>(r) - (80) 

i6 -* 
-g^2 O v -xM+Aix v<£)} 



Ai = A + a: B = rot A; (81) 



is an effective summary vector potential composed from 
the physical external vector potential A generating the 
magnetic field B, and from a gauge vector potential a 
created by the vortices. The effective field A\ represents 
the frustration arising in the system, when the density 
of the particles p(r) — \ip(r)\ 2 deviates away from the 
quantized Lauglin's density n [42] , which determines the 
fractional filling factor v = l/m with m integer. The 
density no is named the flux density being determined 
by the magnetic field B through the magnetic length I in 
the form n = l/m27rP, where I 2 = Hc/eB. The equation 
of motion for vector A\ in a static case is: 



0V x Ai = (ip*ip - n ); = 2ir/m 



(82) 



The proposed phenomenological G-L theory allows us to 
understand that the creation energy of a single vortex is 
finite and that the vortex has a fractional charge. The 
difference between the FQHE and ordinary superfluidity 
was explained by the strong phase fluctuations induced 
by the frustration. Zhang, Hansen and Kivelson [42] de- 
rived their field-theory model starting from the micro- 
scopic Hamiltonian. They constructed the G-L theory in 
a way similar to Girvin but contrary to Girvin in their 
approach the Chern-Simons term contains only the gauge 
field a(r) [42]. As in the previous papers [53] it was con- 
firmed that the disturbances of the localized density mov- 
ing the system away from the good filling fractions lead 
to creation of single-particle excitations. These quasipar- 
ticle and quasihole excitations have the form of vortices 



with static nonuniform finite-energy solutions. Side by 
side with the single-particle excitations in the Ref. [42] the 
collective elementary excitations were discussed. For this 
end the Lagrangian was expanded up to terms quadratic 
in 8<p and 5a about the constant solutions corresponding 
to vacuum expectation values. The fluctuating values 5<fi 
and 5a were represented in the form of plane waves with 
vector q. The dispersion relation was found in the form 
[42] 



uj 2 (q) = (enB) 2 + ^nq 2 {nq 2 



SXno) 



(83) 



It has a gap in the point q = proportional to the 
external magnetic field B. For negativ A, but for suffi- 
ciently small parameter |A|/k the dispersion curve has 
a roton-type behavior with the same shape as was de- 
rived by Girvin, MacDonald and Platzman [54]. The 
GL theory proposed by [42] describes the incompress- 
ibility, fractional charge and fractional statistics of the 
quasiparticles. But being a coarse-grained version of the 
FQHE it makes errors on the magnetic length scale. It 
treats the gauge field with a mean-field approximation, 
and reproduces correctly the long-wavelength effects of 
the quantum Hall systems excluding such details as the 
description of the vortex core. The idea that the long- 
wavelength effects of the physical magnetic field are can- 
celed by the gauge field was also suggested by Laughlin 
[55] and in Ref. [56]. 



VIII. POINT VORTICES UNDER THE 
CONDITIONS OF FQHE 

Because the vortices play an important role in the un- 
derstanding of the FQHE we will provide here more in- 
formation on this subject. The presentation below be- 
ginning with classical hydrodinamics and proceeding to 
the quantum vortices is given following the paper [57] 
and Enger [44] and Myklebust theses [58] . An ideal fluid 
without viscosity is described in classical hydrodynamics 
by the continuity equation. 



dp 
dt 



+ v(H = o 



and Euler's equation 



dv , ^. Vo 
- + (,.v> = -- 



(84) 



(85) 



where p, p and v are the density, pressure and velocity 
field correspondingly in each point of the liquid. The vor- 
ticity is defined in 3D hydrodynamics as Co = V x v. If the 
liquid is not only ideal but also isentropic with constant 
entropy along it, then the vorticity uj obeys to a sup- 
plementary continuity equation. The flow is irrotational 
with a potential flow if <2 = at all points of the fluid. 
In this case one can introduce the velocity potential 4> 



V x v = 



(86) 
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In physical fluids the vorticity is localized in small areas. 
Outside the vortices most of fluid is irrotational. In a 3D 
liquid the vortex is a tube with the strength k defined as 



C3da 



vdl 



(87) 



The Helmholtz theorem (also known as Kelvin's circula- 
tion theorem) says that in the absence of rotational ex- 
ternal forces a fluid that is initially irrotational remains 
irrotational all the time. In case of 2D fluid the notion of 
point vortex with zero area is introduced. The velocity 
field generating such a vortex may be represented by the 
expressions 



2nr w 2tt V r 2 
eg = j cos 9 — i sin 9; e r = i cos I 
_ -0 -0 d _ 
ox Oy or 



-> x 



+ j sin t 
1 d _ 
r oO 



Here n is the vortex strength, whereas the unit vectors 
i, j, e r and eg corresponds to rectangular and polar 2D 
coordinates. Following Ref [41] we must take into account 
the definition of the curl in the 2D space, namely that 
the curl of the vector is a scalar and the curl of the scalar 
is a vector as follows 

oj = Curlv = Vxv = e ij d iV j\ (CurlSf = e lJ d 3 S (89) 

where e lJ is an antisymmetric tensor with the properties 

£ 12 = _ £ 21 = 1; £ 11 = £ 22 = 

These rules lead to the vorticity of the point vortex with 
the velocity field 



(90) 



lo{t) = Curlv — — Alnr = K<5^ 2 '(r) 
2tt 



The velocity field created by a point vortex has a singu- 
larity. It is irrotational or potential almost in all space 
except of the origin in the point r = 0. By this reason the 
vortex area is zero. Nevertheless the summary vorticity 
due the singularity (90) is finite. In the same way the 
circulation of the vortex is also finite as follows 



/u)d 2 r= lim <p -^—egdl — k; dl — rdOeg 
r^oo J 2nr 



(91) 



A fluid containing a point vortex will have potential flow 
almost everywhere. A point vortex in an incompressible 
liquid has energy 



/ 
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err = In — 

2 2tt a 



(92) 



where R is the length scale of the whole system and a 
is the core radius. A classical system of N point vor- 
tices in an incompressible liquid has the kinetic energy 
associated with each vortex and the interaction energy 
between them. This interaction does not come from an 
electric charge of the vortices because they are neutral. 



For two vortices with guiding centers R\ and R2 it is use- 
ful to define a guiding center of a pair R gc and its relative 
coordinate R re i in the form 



R 



R\ + R2; Rrel = i?l — i?2 



(93) 



The equations of motion for a pair of vortices with equal 
strengths K\ = k 2 = k are 



KY r el 



rel 



(94) 



These equations describe a circular motion around a fixed 
point named as stationary guiding center with an angular 
velocity ? depending on the constant separation distance 



of the vortices 



R 



red 



as follows 



7T R re l 



(95) 



For a pair of vortices with opposite vorticities k = Ki = 
— K 2 i.e. for a vortex-antivortex pair the equations of 
motion are 



Xg C — 



K 



Yr, 



R 
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rel 



, X r el ~ ^reZ 



(96) 

The vortices will not move relative each other, but will 
follow a straight line perpendicular to the vector R re i 
connecting the vortices [44]. 

This picture is exactly the same as the structure of a 2D 
magnetoexciton moving with wave vector k perpendicu- 
lar to the vector d connecting the electron and hole in 
the pair with a constant distance d — kl 2 at a given k. 
One can remember that the existence of the quantum 
vortices was suggested for the first time by Onsager [59] , 
who proposed that the circulation in the superfluid He-II 
is quantized with the quantum of circulation h/m. The 
quantum vortices in the He-II were discussed by Feynman 
[60], whereas a quantized line was observed by Vinen [61]. 
The quantization of the vorticity in the He-II can be ex- 
plained in the frame of GL theory. The velocity field of 
a superfluid described by the wave function 



is 



can be written as 



h -a 
v = —S/S 

m 



(97) 



(98) 



The circulation around a close path C becomes 

k = I vdl= — I VSdl = —SS (99) 
/ m J m 

8S is the change in the phase of the wave function, as one 
moves around the close path C. But the wave function 
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must be single valued. By this reason SS must be integer 
multiple of 27r. It means that 



h 

k = — 2tts, s = 0,±1,±2... 

m 



(100) 



The vorticity of the quantum vortex has discrete values 
with the quanum h/m. This definition of vorticity differs 
from the classical hydrodynamics co = V x v. The only 
rotational invariant wave function having the property 
(100) being written in polar coordinated has S = s9 



<P{r) = f(r)e 



isO 



(101) 



It produces the same velocity field as the classical point 
vortex 



hs _ k 

v = ee = -z — e e 

mr 2irr 



(102) 



The kinetic energy 

E = J^mv 2 d 2 r= ^s 2 \n(R/$ (103) 

is now expressed through the coherence length £ instead 
of the core radius a. The cutoff at £ is used to avoid the 
logarithmic divergence neared the vortex core. 
Inserting the vortex function (101) into the Ginzburg- 
Pitaevskii-Gross equation, Myklebust [44, 58] found the 
following equation for the function f(r) 



dr 2 



10/ 
r dr 



+ 2 



/ - 2/ 3 = 0. 



(104) 



It depends only on s 2 . Contrary to the He-II, the Bosc- 
Einstcin condensate in superconductors is formed by the 
Cooper pairs with the charge q = — 2e instead of the neu- 
tral atoms. The type II superconductors allow the mag- 
netic field to penetrate in metals forming quantized vor- 
tices, while in the type-I superconductors the magnetic 
field cannot penetrate [62]. The quantized vortices exist 
in the form of filaments named Abrikosov's lines. They 
have a mixed electron and electromagnetic field origin 
and were described for the first time by Abrikosov [63] on 
the base of the G-L theory with nonzero electromagnetic 
field A. It was shown that the magnetic flux through the 
vortex tube is quantized with the flux quantum tfio 
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(105) 



The total energy per a unit length of the vortex tube is 
finite and equals to 



(106) 



where A is the penetration length of the magnetic field 
into the II-type superconductors as was introduced by 
F. London and H. London [64] and £ is the correlation 



length between the electrons in the Cooper pair. Girvin 
[52] suggested that contribution of the electromagnetic 
field in the resultant current density j(r) determined in 
the case of FQHE as 

3(f) = \ {^*(r)(-i/lV)V(r) + V(r)(i7iV)V*(r)} 

+ -A^*(rU(r) (107) 
c 

reorganizes the point vortex state in such a way that its 
resultant circulation at great distance r — > oo will be zero 



j(r)dl = 



(108) 



It is possibly only for the condition when the magnetic 
flux through the vortex surface is quantized in the form 



J rotAd 2 f = <j> AdT= -r 



2i:hc 

vol A<l 2 r= (J) Adl = -m0 o ; <f>o = -j^p (109) 



This value being multiplied by n \e\ jc compensates ex- 
actly the circulation arising from the electron part of the 
current density 

j> \ {*P*{r)(-ihV)iP(r) + V(r)(i/iV)V*(r)} dl 

= 27rft™n (110) 
because the wave function ip(r) has the form 

V'M = V^f(r)e lm8 ; f(r) -> 1; r -> oo (111) 



The number of magnetic flux quanta — m must be op- 
posite to the magnetic quantum number of the electron 
wave function. The creation energy of such point vortex 
is finite and no extensive as in the case of a pure electron 
vortex. As was mentioned by Girvin and MavDonald [53] 
the isolated vortex cost only a finite energy. They can be 
excited thermally by one. Earlier it was necessary to cre- 
ate a pair vortex-antivortex with finite creation energy 
for a pair as a whole, but with extensive energy for each 
of them. Only in the last case the Kosterlitz-Thoulless 
phase transition was possibly being related with the un- 
binding of the vortices in the pairs. 



IX. GAUGE TRANSFORMATIONS AND 
STATISTICAL GAUGE FIELD 

Girvin and MacDonald [53] revealed a hidden symme- 
try of the Laughlin's [65] ground state wave function de- 
scribing the FQHE of the 2D one-component electron gas 
(OCEG). This wave function is 



ip(zi,...,z N ) = Y\ (zi - Zj) m exp 



i<3 



Zk\ 



(112) 



The filling factor of the lowest Landau level (LLL) is 
a fractional integer v = 1/m, with integer m > 1. 
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%k = x k + ij/feare the complex coordinates of the par- 
ticles in symmetric gauge. With respect to the inter- 
changing of any two particles the wave function (112) is 
anti-symmetric at odd values of m and symmetric at even 
values, describing the fermions and bosons, correspond- 
ingly. But changing the phase of the wave function (111) 
using a singular gauge transformation 



1pnew(zi, -,Z N ) = CXp 



X7p(zi,...,Z N ) =Y[\ 

i<j 



■ ^2 ar g ( z i 

i<3 



Zj) 



cxp 



-J]Tm 2 ( 113 ) 



we have obtained a bosonic type wave function at any in- 
teger values of m > 1. The off-diagonal matrix elements 
of the density matrix p(z, z') calculated with the function 
(111) arc short-ranged with a characteristic scale given 
by the magnetic length, while those calculated with the 
wave function (112) p(z, z') have a slowly decreasing be- 



havior with a power law | z 



</2 



. The singular gauge 



density matrix p{z, z') has an off-diagonal long-range or- 
der (ODLRO). The physical origin of this difference is 
related to the presence of the vortices induced around 
each particle under the influence of the magnetic flux 
quanta, as was explained by Stormer [66] . 

The presence of the vortices can be demonstrated using 
more simple example proposed by Enger [44] with a wave 
function ij){z) of two particles depending only on their 
relative coordinate z. It is supposed that ip(z) obeys 
to any statistics and after the particle interchanging it 
becomes 



ip(e i7r z) = e l0 iP(z); 



9 = 7r(2n + 1) for fermions 
9 = 2irn for bosons 



0,±1,±2... 



(114) 



A gauge transformation e lr> ^ transforms the wave func- 
tion ip{z) into another bosonic type wave function 

e lv{z) ^(z) = 4>(z) = 4>{e™z) (115) 
To satisfy this requirement and the equalities 

e"'( e " rz )v> (e t7T z) = e l1liz mz) = e^ e " z ) e w ^{z) (116) 
the function rj{z) must satisfy the equation 

9 9 y 
7](z) — r\ (e t7T z) = 9; r](z) = arg z = arctan — 

7T 7T X 

(117) 

The transformation of the wave function (115) must be 
accompanied by the transformation of the electromag- 
netic field potentials [44] 



i A ^i A » +d ^ z) = 

= -^{A^ + a^; p = 0,1,2 



(118) 



In such a way side by side with the electromagnetic po- 
tential Afj_ supplementary gauge potential a M created by 
the vortices appears. 



he 



aJr) = —dnrjtr) = -— <9u arctan- (119) 

e 7re x 



The statistical gauge vector potential has the expression 



a = Curl In r = V x Inr 

7re 7re 

a-i = ^^-e^dj lnr; i, j = 1,2 
7re 



(120) 



This vector potential is created by the vortices arising 
near each particle. It leads to magnetic field strength 
[41] 

b(r) = Curla(r) = V x a(r) = e^didj = 

= A In r = 5 (2) r); (121 

en e 
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= <5( 2 )(r); 



The magnetic flux created by this magnetic field is 



/ 
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b{r)d 2 r = = O ; 0o = — (122) 

e 7r e 



It equals to — (2n + l)<j>o, when the initial particles de- 
scribed by the function tp(z) were fermions, and equals to 
—2n(j)Q for boson wave function ip{z). This result shows 
that the initial fermion particles each of them attaching 
an odd number of flux quanta transform themselves into 
a composite bosons described by the new wave function 
<t>(z) which obeys to Bose statistics. The effective mass 
m and charge e remain the same at least in the given 
approximation but the composition and statistics of the 
final quasiparticles are changed. It is said that the elec- 
tron attached an odd number of flux tubes, though in 
reality such tubes do not exist. We can say that in our 
case the initial particles are fermions or electrons whereas 
the final quasiparticles are bosons. Formula (121) may 
be generalized for any number of particles, which create 
in a common way the resultant magnetic field 

,. . 2hc6 v— a ,.9 . _ 28hc , . 

b(r) = y« 2 (f-f,) = p(r) (123) 



where p{r) is the density of the particles. 
As was mentioned above, Zhang, Hanson and Kielson [42] 
have generalized the Ginzburg-Landau theory introduc- 
ing into the Lagrangian a supplementary Chern - Simoms 
[35] term related with the influence of the statistical 
gauge field. The Lagrangian of the Ginzburg-Landau- 
Chern-Simons (GLCS) theory in the form presented by 
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Enger [44] looks as 



L = ihcj)* [d t + — (A + a ))(j) + 



(124) 



Here the following denotations were used: p, v 1 a = 
0,1,2; do = d t , di = {di = d x , 82 = d y }. The tensor 
£ p.ucT nag nonzero components only for different values 
of p,u,a. They change the signs at any permutations of 
two indexes as follows: 



e 012 = 1, e 021 



-1, 



(125) 



.102 



-1, e 120 = 1 etc. 



The external electromagnetic 2D vector potential A and 

the scalar potential A are taken as A^ = (^A a ,Aj. 

is the statistical gauge potential with three components. 
Two of them a = (01,02) generate the statistical "mag- 
netic" held and the third component ao gives rise to the 
statistical "electric" field. Two parameters m and e of 
the Lagragian (124) are the effective mass and charge of 
the final type quasiparticles obeying the Bose statistics. 
They can differ from the mass and charge of the initial 
particles. The A and po parameters are typical for the 
G-L theory, while p is the Chern-Simons parameter [67]. 
Variations of (51) with respect to <fi* leads to the nonlin- 
ear Shrodinger equation: 



ihd t (A + a ) 

c 



~ 2m 



V + 



(126) 



-A(^-poH 
The variation of (49) with respect to a gives 

pe v didj — e<j)*4> = ep 

which can be transcribed as 

pCurla = pV x a = pb = ep 



(127) 



(128) 



Parameter p can be determined comparing it with the 
expression (119) 



26hc 



p 



(129) 



For the initial fermion particles with 9 = 7r(2n+l) the pa- 

2ithc 

ramctcr p of the Lagrangian (124) equals to ^— (2n + 

e A 

1) = — (2n+ 1) — . The third variation with respect to 

e 

di gives 



pe lJ ( Va - d t a) = ej(r) 



(130) 



where j(r) is the current density 2D vector determined 
as 

H 



v + ^(l+a))<rj (LTD 



The equation (130) states that the statistical "electric" 
field 



1 d _ 

e = -Va + -—a 
c dt 



(132) 



is related with the particle current density j(r). 

The energy density of the GLCS system in the purely 

static external magnetic field (A = 0) equals to 
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2m 
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A simple solution to these equations can be obtained by 
setting 



9( r ) = VPo e 
It must obey the equation 



iS(r) 



In the case S = Const we have 
A + a = 



(134) 



(135) 



(136) 



It means that the corresponding magnetic fields B = 
V x A and fc=Vxa cancel each other. The final quasi- 
particles named as composite bosons [68-70] do not feel 
the net magnetic field and behave as bosons interacting 
with each-other via the <fi 4 type interaction. The notion of 
composite particles (CPs) consisting from electrons and 
attached magnetic flux quanta was introduced firstly by 
Wilczek [68]. 

As was mentioned above, the collective elementary exci- 
tations of the described ground state are the plane waves. 
Their dispersion law has a gap, and this means that the 
system is an incompressible quantum liquid, which can 
not be excited by a very small perturbation [69, 70]. 

Above we have discussed the case when the initial wave 
function ip(z) with Fermi statistics was transformed into 
another wave function <p(z) obeying Bose statistics us- 
ing a singular gauge transformation. Read [40, 71-73] 
investigated the system of 2D charged bosons interact- 
ing with a transverse magnetic field and between them- 
selves. The filling factor of the LLL was supposed to be 
one. It means that there is one flux quantum for each 
particle. Following Read [40, 71-73] it is equivalently to 
say that there exists one vortex for each particle. In this 
case the vortex has the charge of opposite sign in com- 
parison with the charged boson and the Fermi statistics. 
Now the gauge transformation attaching one vortex to 
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each charged boson will create composite particles with 
resulting charge zero and with Fermi statistics. The neu- 
tral composite fcrmions (CFs) will move in zero mag- 
netic field. Such system can be described in the frame 
of the Fermi-liquid-theory. Another variant proposed by 
Halperin, Lee and Read [39] was considered in Section 5. 
The starting Hamiltonian describes the electrons form- 
ing a 2D electron gas (2DEG) with filling factor v of the 
LLL equal to one half (y = 1/2). Now for each elec- 
tron there are two flux quanta or two vortices each of 
them having the charge — e/2 and Fermi statistics. Two 
vortices are equivalent to one 2-fold vortex with charge 
— e and Bose statistics. The gauge transformation of the 
wave function will transform the initial charged electrons 
into the composite neutral fermions each of them consist- 
ing from one electron and 2-fold vortex. The Hamilto- 
nian of the system will be changed because instead of 
the external magnetic field will appear a supplementary 
gauge magnetic field, which in well definite conditions 
cancels exactly the external magnetic field. The initial 
charged fermions were converted in neutral CFs moving 
in a zero resulting magnetic field. The fictitious Chern- 
Simons "magnetic" field created by the vortices being 
averaged in the mean-field approximation cancels exactly 
the external magnetic field only in statistical sense and 
under the definite conditions. It happens when the mean 
density of the electrons corresponds to the fractional in- 
teger filling factor. In the present example with v = 1/2 
the gauge transformation does not modify the statistics 
of the composite particles (CPs). As earlier, they are 
neutral CFs in a zero magnetic field. The singular gauge 
transformations were firstly introduced by Wilczek [68] . 

The single-particle elementary excitations appear in 
the form of the fractionally charged vortices. They are 
fcrmions and have finite creation energy as was under- 
lined by Girvin [52], and Girvin and MacDonald [53]. 
Read [39, 40, 71-73] argued that the ground states of the 
systems in the condition of FQHE with different frac- 
tional integer filling factors v = 1/m with m = 1, 2,3, ... 
contain electrons bound to vortices, since such binding 
lowers the system's energy. A m-fold vortex carriers a 
charge —evm in the fluid, where e is the electron charge 
e = — |e|. The electron- m- vortex composite, named as 
CP, at v = 1/m has a net charge zero and behaves like 
a particle in a zero magnetic field. The vortex is sensi- 
tive to the density of electrons, which can vary in space 
and time even when the external magnetic field and the 
average filling factor are fixed. The m-fold vortices are 
fermions for m odd and bosons for m even. The compos- 
ite boson particles can undergo the Bose-Einstein con- 
densation (BEC), because it minimizes their "kinetic" 
energy. Just the BEC of CBs is the interpretation of the 
Laughlin's states [65] . The origin of the "kinetic" energy 
is the potential energy of the interaction between the par- 
ticles. In the case of electrons it is the Coulomb electron- 
electron interaction which is not canceled by the gauge 
transformation and CS gauge potential. It is named as 
"kinetic" because it depends on the wave vector of the 



operators (139). The bound objects such as CPs do, in 
fact, have such an effective "kinetic" energy. There is 
an attraction between an electron and m-fold vortex. It 
plays for the electron the role of a correlation quasihole. 
As was shown in (139) the CPs may exist in the form of 
plane waves and the many-particle wave functions also 
can be characterized by the wave vector k. 

The creation operator in the coordinate representation 
(56) can be rewritten in momentum representation as 
follows 



V> T (f) 



(137) 



A CP with k — would have the electron exactly at the 
zeroes of the wave function or in the center of the vortex, 
whereas the CP with wave vector \k\ ^ has the electron 
displaced by the distance |fc| I 2 from their center. One 
can say that the electron and its correlation quasihole 
or in another words the electron and the m-fold vortex 
experience a potential V (|fc|) due to the Coulomb inter- 
action of the electron with other electrons excluded from 
the vortex core. All these interactions take place in the 
presence of the neutralizing background. The electron 
and the m-fold vortex experience the magnetic field of 
the same strength. Both components of the pair drift in 
the same direction perpendicular to the vector connecting 
their centers, so that their separation remains constant 
and equal to \k\l 2 . The energy of a pair is V (\k\) and its 
group velocity is dV (\k\) /d\k\ [40, 71-73]. Wc can add 
that this picture coincides with the structure of the 2D 
magnetoexciton, where the energy V(|fc|) equals to the 
expression E (| k |) [39] 



E{\k\) = 2^W Q Sin 2 
Q 
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(138) 



Here eo is the dielectric constant and S is the layer surface 
area. 

Girvin, MacDonald and Platzman [54] elaborated the 
theory of the collective elementary excitation spectrum 
in the case of FQHE, which is closely analogous to the 
Feynman's theory of superfluid helium. The predicted 
spectrum has a gap at k — and a deep magneto-roton 
minimum at finite wavevector, which is a precursor to the 
gap collapse associated with Wigner crystal instability. 
They supposed the existence of only one branch of the 
collective elementary excitations spectrum. In this ap- 
proximation named as single mode approximation (SMA) 
they have constructed the wave functions of the excited 
states 4*k acting with the operator of the particle den- 
sity pk on the ground state wave function ip g in the form 
4>k = Pkipg- They determined the energy of the excited 
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state A(k) as 



A(fc) 



t \(H - Eo)\4> k ) 



(4>k\<i>k) 



s(k) 



(139) 



where f(k) is the oscillator strength and s(k) is the static 
structure factor. The total oscillator strength sum is sat- 
urated by the cyclotron mode contribution, and f(k) has 
a dependence of the type \k\ 4 . As was established by Lee 
and Zhang [74] the influence and the contribution of the 
quantum vortices to the dynamical and static structure 
factors is important. It leads to dependence s(k) <~ \k\ 4 
at k — > 0. In this case A(fc) has a gap. Neglecting the 
influence of the quantum vortices the dependence s(k) is 
proportional to |fc| 2 and the energy spectrum is gapless 
A(fc) w k 2 at k — > as a Goldstone mode [2]. In conclu- 
sion, for the FQHE in the case v < 1 with fractionally 
filled Landau level the Pauli principle no longer excludes 
the low-lying intra-Landau-level excitations. They exist 
side by side with the inter-Landau-lcvel excitations. The 
last excitations have a cyclotron energy gap [75]. 

We are studying a coplanar electron-hole (e-h) system 
with electrons in conduction band and with holes in va- 
lence band in a strong perpendicular magnetic field. Pre- 
viously such system has been studied in a series of papers 
[9, 22-25, 27, 76]. Most of them were dedicated to the 
theory of 2D magnetoexcitons. All the same, there were 
papers dedicate to another aspects of these systems. For 
example MacDonald, Rezayi and Keller [77] as well as 
Joglckar and MacDonald [78] have discussed the photo- 
luminescence (PL) spectrum in the FQHE regime. It 
was mentioned that the PL spectrum does not exhibit 
anomalies associated with the FQHE. However when the 
electron and hole layers were separated a new peak in 
the PL spectrum appears, when the filling factor exceeds 
a fraction vq at which an incompressible quantum liquid 
occurs. The new peak is separated from the main spectral 
features by the quasiparticlc-quasihole gap. We are inter- 
ested in the distribution of the flux quanta in the case of 
e-h system with equal average numbers of electrons and 
holes N e = Nh with the filling factor v = N e /N, where 
TV is the total number of flux quanta , where S is the 
layer surface area and 27r/ 2 is the area of the cyclotron 
orbit. For the fractional integer filling factor there are an 
integer number of flux quanta per each e-h pair. The cre- 
ation of the vortices in this case is not studied at present 
time. But one can expect that in the case of magnetoex- 
citons they will be neutral, whereas in the case of pure 
electron and pure hole vortices their "magnetic" gauge 
fields will compensate each other, so that the charge, the 
statistics of the particles and the external magnetic field 
will remain the same in the mean-field approximation 
with equal densities of electrons and holes. Neverthe- 
less due to quantum fluctuations and the deviations in 
space and time of the electron and hole densities from 



their average values one can expect the influence of the 
pure electron and hole quantum vortices on the physics 
of magnetoexcitons side by side with the influence of the 
neutral quantum vortices formed by the magnetoexcitons 
themselves. The last quantum vortices determine the 
Berezinskii-Kosterlitz-Thouless phase transition [79, 80]. 



QUANTUM HALL EXCITONS IN BILAYER 
ELECTRON SYSTEMS 



In this section we give a short review of the Bosc- 
Einstein condensation (BEC) of the quantum Hall exci- 
tons(QHExs) arising in the bilayer electron systems un- 
der the conditions of the quantum Hall effect (QHE) at 
one half filling factor v = 1/2 for each layer and the total 
filling factor for two layers equal to unity u t = 1. The 
purpose is to compare this phenomenon with the case 
of BEC of two-dimensional (2D) magnetoexcitons. Such 
comparison will give a better understanding of the un- 
derlying physics and allows to verify the accuracy of the 
made approximations. 

In the Ref.[81] Fertig investigated the energy spectrum 
of a bilayer electron systems in a strong perpendicular 
magnetic field and introduced the concept of the inter- 
layer phase coherence of the electron states in two ad- 
jacent layers, which leads to the model of quantum Hall 
excitons under the condition of their BEC. Unexpectedly 
a strong evidence of exciton BEC was ultimately found 
in such surprising system as a double layer 2D electron 
system at a high magnetic field [82] . In the QHE regime 
the excitons consist of electrons in the lowest Landau 
level (LLL) of the conduction band of one layer being 
bound to the holes which appear in the LLL of the con- 
duction band in another layer. The formation of such 
unusual holes is due to the possibility to consider the 
half-filled LLL by electrons of the conduction band, for 
example, of the first layer as being completely filled by 
electrons with filling factor v = 1 and simultaneously 
being half-filling by holes in the same conduction band. 
The full-filling electrons of the first layer are considered 
as being compensated by the impurity doped adjacent 
layer and the theoretical model takes into account only 
the holes in the first layer and the electrons in the sec- 
ond layer. Both components belong to the LLLs of the 
same conduction band and are characterized by a half- 
filling factor for each of them. This new type of excitons 
named QHExs appears whenever the temperature and 
the layer separation are small enough and the total den- 
sity rit of electrons in the double layer system equals to 

, i eB 1 

the degeneracy - — = n t = — - 

he ' 



where I is the magnetic 

length. The total filling factor v t — n t 2Trl 2 equals unity. 
The new collective electronic state introduced by Fertig 
[81] exhibits several dramatic electrical transport prop- 
erties revealed in Ref. [83-85]. As was mentioned in [82] 
the BEC of the QHExs reflects the spontaneously broken 
U(l ) symmetry in which the electrons are no- longer con- 
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fined to one layer or to the other, but instead they reside 
in a coherent linear combination of the two layers. This 
interlayer phase coherence develops only when the effec- 
tive interlayer separation d/l is less than a critical value 
(d/l) c . At large d/l the bilayer system behaves qualita- 
tively like the independent 2D electron systems. Follow- 
ing [86] this new state can be distinguished as a Fermi 
liquid state of composite fermions. It is unique because 
unlike other QH states it possesses a broken symmetry in 
the absence of the interlayer tunneling. It can be viewed 
as a pseudospin ferromagnet with the pseudospin encod- 
ing the layer degrees of freedom or as an exciton BEC 
with QHExs formed from electrons and holes confined to 
different layers. 

There are two energy scales in the double-layer systems. 
One is the potential energy V between the electrons in 
different layers. The second is the energy gap A$as be- 
tween the symmetric and asymmetric states of electrons 
in two layers measuring the tunneling amplitude between 
them. The capability of tuning the strength of the inter- 
layer interaction by changing the gate voltage provides 
the opportunity to explore the v t = 1 system through 
its transformation between the weak and strong inter- 
action limits and to study the phase transitions between 
the compressible Fermi liquid and the incompressible QH 
states as a function of d/l [86]. In the most theoretical 
investigations of the QHExs, except the paper [86], the 
simplifying assumption of the fully spin polarized elec- 
trons was used. Below, in our discussions the Zeeman 
energy will be not included. Following the Ref. [87], in 
the absence of the interlayer tunneling there are two U(l) 
symmetries. One is associated with the conservation of 
the total electric charge N\ + N 2 , where Ni and N 2 are 
the numbers of electrons in two layers, and the other is 
related with the conservation law of Ni — N 2 . For these 
conditions the gapless mode appears. It is the Nambu- 
Goldstonc(NG) mode arising from the broken U(l) sym- 
metry associated with N\ — N 2 and characterized by the 
off-diagonal long-range order in the tunneling operator 
a\ p a 2q , where a\ p and a 2q are the electron annihilation 
operators in the LLLs of the conduction band of two lay- 
ers. 

Within the mean-field effective theory the appearance of 
the gapless mode may be attributed to the coherent fluc- 
tuation of the electron flux and density describing the 
relative fluctuations of the electron densities in two lay- 
ers. At finite interlayer tunneling the number Ni — N 2 is 
no-longer conserved. As well, the currents in each layer 
are no-longer separately conserved. 

The collective excitation spectrum of the two-layer elec- 
tron system with v t — 1 was investigated by Fertig [81] 
on the basis of the theoretical model without tunneling 
but with different interlayer separation including d = 
and taking into account that at d > the Coulomb in- 
terlayer electron-electron interaction is smaller than the 
intralaycr interaction. The ground state wave function 
proposed by Fertig [81] introduces the interlayer phase 
coherence reflecting a new state, in which the electrons 



are no-longer confined to one layer or to another, but in- 
stead they reside in coherent linear combinations of the 
two layer states as follows 



n { ua u + v 4t) io) , 



u 2 + v 2 = 1 



(140) 



The lowest levels of the Landau quantization in the Lan- 
dau gauge are characterized by the quantum number 
n = and the uni-dimensional wave number t, with 0) 



being the vacuum state. The equality u 2 



1/2 re- 



flects the half-filling of the LLL in each layer. Introducing 
the hole operator d\ , d t for the first layer instead of the 
operators a\ t and a\ t the function (144) was transcribed 
in the form 

M =U{ u + va t d -t) M>>. 
t 

t 

o-2t = a, t , a 2t — a\, ait — d_ t , a lt = d-t 

The operators a\d}_ t create the electron-hole pairs with 
total wave vector equal to zero. The wave function (141) 
can be interpreted as describing the BEC of the QHExs. 
This model is similar to the case of BEC of 2D magnc- 
toexcitons studied in Refs. [22-24, 76]. In the last case 
the holes were formed in the frame of the valence band. 

The valley-density two-particle integral operators in- 
troduced in Ref. [81] in the electron-hole representation 
are 



= Zy«»"y qx dq x +a\ x aq x ± 
±(dq x aq x -a\ x d\ x )]; (142) 

y~ ( y +t -Y +t -y-' 

PM =y^e iq v tl2 (dq x aq x 



t 



t ^\x d \x ) 

2 T +t -y + * -y- ( 



p F (g) = ]T e ^V 
t 



J „ dq x 



e*"" la q x aq x - u Qx 

y +t y + * -y-' y 



We introduce our designations for the exciton, optical 
and acoustical plasmon operators with holes in the con- 
duction band, as follows 
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In the case of holes in the valence band there are opposite 
signs in the expressions for p(q) and D(q). The relations 
between two sets of operators are 



--(q) = p(q)±VN(d(q)-dH-q)) 



•?)) 



M) = VN(d(q)+dH- 
p F (q) = D(q) 



(144) 



The response functions were introduced as follows [81] 



X±(«,w) = -i / dte iut ([p^(q,t),p ± (-q,0)]) 



Xz(q,uj) = -iJ dte^ (\p z (q,t),p z (-q,0)}) (145) 
o 

oo 

XF(q,u>) = -i J dte^ 1 ([p F (q,t), p F (-q,0)]) 
o 

The poles of these functions represent the excitations of 
the system. The excitations may be thought as a valley- 
density waves or pseudospin density waves. The calcu- 
lation of the response functions were effectuated by Fer- 
tig [81] using the diagrammatic expansion elaborated by 
Kallin and Halperin [75]. The response functions were 
written in terms of vertex functions, Green's functions 
and self-energy parts. This approximation is shown by 
the diagrams in Fig 3a Ref. [81] neglecting the diagrams 
which contain supplementary the c-h bubbles. Their con- 
tribution is negligible only when the excited Landau lev- 
els (ELLs) are taken into account and the bubbles have 
an energy denominator huj Cl where uj c is the cyclotron 
frequency increasing in the strong field limit. A self- 
consistent calculation of the vertex function including the 
bubbles is quite difficult. Below we will present the re- 
sults obtained in Ref. [81] for the energy spectrum of the 
collective elementary excitations. 

For d = the interaction Hamiltonian is invariant 
under the unitary transformation SU(2). The specific 
choice of the ground state wave function (141) is a bro- 
ken symmetry state and one expects the appearance of a 
NG mode. At d = the NG mode has a dispersion re- 
lation uj(k) ~ k 2 for the small k. For d > the problem 
can be mapped onto an equivalent spin system with linear 
dispersion relation at small wave vectors. The NG mode 
at d > has a linear dispersion law with a slope depen- 
dence on d, which is similar with that of the acoustical 
plasmon mode of a two-layer system in the absence of the 
magnetic field [88]. To better understand this result one 
may recall the BEC interpretation of the ground state 
wave function. Indeed, at d > the inter-layer electron- 
hole Coulomb attraction is smaller than the intralayer 
electron-electron and hole-hole repulsions, that leads to 
a resultant repulsion in the system and to the transfor- 
mation of the parabolic dispersion law into the linear 
one at small values of wave vectors as in the Bogoliubov 



theory of superfluidity of the Bose gas [1]. At kl of the 
order unity the dispersion law of Ref. [81] develops a dip 
at certain critical values of d, indicating that the system 
trends to undergo a phase transition. Another consider- 
ations concerning of the gapless modes in the FQHE of 
multicomponent fermions can be found in Ref. [89, 90]. 
The above branch of the energy spectrum corresponds 
to the response function x±(<li UJ ). The operators p ± (q) 
describe two superpositions of the optical plasmon and 
exciton mode operators. There are two other operators 
p z {q) and Pf(q) which describe the pure exciton modes 
and the acoustical plasmon mode. As was established 
in Ref. [81] the last two response functions x z (q,cj) and 
Xf(<7>w) have no poles when the excitations in the frame 
of the LLLs are considered. The excitations associated 
with these functions are considered to be higher in en- 
ergy than the NG mode discussed above by an amount of 
energy of the order haj c . It means that the pure exciton 
and acoustical plasmon modes in the system of BEC-ed 
QHExs cannot be described by the NG gapless modes. 
As will be shown in the next section, in similar case of the 
BEC of coplanar 2D magnetoexcitons the optical plas- 
mon branch is also the unique NG mode. The exciton 
branches (energy and quasienergy) of the spectrum have 
the gaps in the point k = 0, the roton-type behavior at 
intermediary values of the wave vectors and a saturation 
dependences at k — > oo. At the same time the acous- 
tical plasmon branch in the case of magnetoexcitons in 
the range of small wave vectors reveals the absolute in- 
stability. Its values are pure imaginary. In the case of 
BEC of magnetoexcitons there is one NG optical plas- 
mon mode, two gapped exciton modes and one unstable 
acoustical plasmon mode. This agrees qualitatively with 
the results obtained by Fertig [81] in the case of BEC of 
QHExs where a single gapless NG mode of the optical 
plasmon type was revealed while the other modes of the 
spectrum were not identified at infinitesimal energies. 



XI. TRUE, QUASI AND UNSTABLE 
NAMBU-GOLDSTONE MODES OF THE 
BOSE-EINSTEIN CONDENSED COPLANAR 
MAGNETOEXCITONS 

In this section we will present the results following the 
Refs. [30, 31] for the energy spectrum of the collective 
elementary excitations arising above the ground state of 
the Bose-Einstein condensed coplanar magnetoexcitons. 

The full Hamiltonian describing the interaction of elec- 
trons and holes lying on the LLLs is: 



H — Hcoui + Hsuppi 



(146) 



Where Hcoui is the Hamiltonian of the Coulomb inter- 
action of the electrons and holes lying on their LLLs: 

Hcoui - W Q \p(Q)p(-Q) -N e - N h 



-p e N e - p h N h (U7) 
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and Hsuppi is the supplementary indirect interactions be- 
tween electrons and holes, which appear due to the simul- 
taneous virtual quantum transitions of two particles from 
the LLLs to excited Landau levels (ELLs) and their re- 
turn back during the Coulomb scattering processes. The 
expression for this interaction was obtained in Ref. [25] 
and has the form: 

tf-ppi = \Bi-iN -L £ v(Q)p(Q)p(-Q) 

Q 

±Y U (Q)^QW-Q) (148) 



Q 



Here p(Q) are the density fluctuation operators expressed 
through the electron p e (Q) and hole ph(Q) density oper- 
ators as follows: 



p e (^) = ]T e ^ 2 a 



Qx a . . 

t+— t-— 

p(Q) = Pe(Q) - p-h(-Q); 
D(Q) = p e (Q) + ph(-Q); 
N e = p e (0); N h = p h {Q)- 



(149) 



N = N e + N h ; W$ 
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The density operators are integral two-particle operators. 
They are expressed through the single-particle creation 
and annihilation operators a p , a p for electrons and fet , b p 
for holes. Here, £o is the dielectric constant of the back- 
ground; p e and ph are chemical potentials for electrons 
and holes, and coefficients V(Q), U(Q) and Bi-i were 
calculated in [25, 29]. 

The starting Hamiltonian (146) has two continuous 
symmetries. One is the gauge global symmetry J7(l) and 
another one is the rotational symmetry 50(2), so that 
the total symmetry is U(l) x S0(2). The gauge sym- 
metry is generated by the operator N of the full particle 
number, when it commutes with the Hamiltonian. It 
means that the Hamiltonian is invariant under the uni- 
tary transformation U{<p) as follows 



U{v)HU-\<p) 



e tfitfi ; [H,N}=0 



(150) 



The operator N is referred to as the symmetry gener- 
ator. The rotational symmetry SO (2) is generated by 
the rotation operator C z (<p) which rotates the in-plane 
wave vectors Q on the arbitrary angle ip around z axis, 
which is perpendicular to the layer plane and is parallel 
to the external magnetic field. Coefficients W^, U(Q) 

and V(Q) in formulas (6) and (9) of Ref. [30] depend 



on the square wave vector Q which is invariant under 
the rotations C z (ip). This fact determines the symmetry 
50(2) of the Hamiltonian (146). The gauge symmetry 
of Hamiltonian (146) after the phase transition to the 
Bose-Einstein condensation (BEC) state is broken as it 
follows from expression (16) of Ref. [30]. In terms of the 
Bogoliubov theory of quasiaverages, it contains a supple- 
mentary term proportional to fj. The gauge symmetry is 
broken because this term does not commute with oper- 
ator N. Moreover, this term is not invariant under the 
rotations C z {<p) , because the in-plane wave vector k of the 
BEC is transformed into another wave vector rotated by 
the angle ip in comparison with the initial position. The 
second continuous symmetry is also broken. Thus, the 
installation of the Bose-Einstein condensation state with 
arbitrary in-plane wave vector k leads to the spontaneous 
breaking of both continuous symmetries. 

We will consider a more general case of k ^ taking 
the case k = as a limit k — >• of the cases with small 
values kl«l. One should keep in mind that the sup- 
plementary terms in Hamiltonian (146) describing influ- 
ence of the ELLs are actual in the range of small value 
kl < 0.5. Above we established that the number of the 
broken generators (BGs) denoted as Nbg equals to two 
{N BG = 2). 

As discussed in previous papers [22-25, 27, 76, 91, 92], 
the breaking of the gauge symmetry of the Hamiltonian 
(146) can be achieved using the Keldysh-Kozlov-Kopaev 
[93] method with the unitary transformation 



d{k))], 



(151) 



where d)(k) and d(k) are the creation and annihilation 
operators of the magnetoexcitons. In the electron-hole 
representation they are [22-25, 27, 76, 91, 92]: 



x J- X 



(152) 



BEC of the magnetoexcitons leads to the formation of 
a coherent macroscopic state as a ground state of the 
system with wave function 



|0)=0 (153) 



Here |0) is the vacuum state for electrons and holes. In 
spite of the fact that we kept arbitrary value of k, nev- 
ertheless our main goal is the BEC with k = and we 
will consider the interval 0.5 > kl > 0. The function 
(153) will be used to calculate the averages values of 

the type (d(Q)D(-Q)^. The transformed Hamiltonian 
(146) looks like: 

(154) 
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and is succeeded, as usual, by the Bogoliubov u-v trans- 
formations of the single-particle Fermi operators 

= ua p - v(p- y)6j^_ p ; 

Op\ip g (k)) =0; (155) 

p p = b (yiC) b P fr (yrC) = 

& - 0; 

Instead of this traditional way of transforming the ex- 
pressions of the starting Hamiltonian (146) and of the 
integral two-particle operators (149) and (152), we will 
use the method proposed by Bogoliubov in his theory of 
quasiaveragcs [1, 51], remaining in the framework of the 
original operators. The new variant is completely equiv- 
alent to the previous one, and both of them can be used 
in different stages of the calculations. For example, the 
average values can be calculated using the wave function 
(153) and u-v transformations (155), whereas the equa- 
tions of motion for the integral two-particle operators can 
be simply written in the starting representation. 

The Hamiltonian (146) with the broken gauge symme- 
try in the lowest approximation has the form 

K = \Y. W Q [P(Q)P(-Q) ~ Me ~ N h - 

Q 

-p e N e - p h N h + ^Bi_iN- 
-^£n<?)/3(Q)p(-Q)- (156) 



BEC of magnetoexcitons with k = are 



-^£tf(Q)£(3)£(-3)- 



Q 



-ijVN (d\k)+d(k)) 

For simplicity another smaller term of this type propor- 
tional to rj was dropped. Here parameter 77, which deter- 
mines the breaking of the gauge symmetry, depends on 
the chemical potential p and on the square root of the 
density, similar to the case of weakly non-ideal Bose-gas 
considered by Bogoliubov [1, 51]. In our case the density 
is proportional to the filling factor v — v 2 and we have: 



\i = n e + n h ; fi = fj, + If, N ex = v 2 N; 
E ex (k) = -It - A(k) + E(k); 
77 = (E ex (k) - n)v = (E(k) - A(k) - fi)v; 

'[KxQ]J 2 ' 



(157) 



E(k) 



2^W Q Sm 2 
Q 



The equations of motion for the integral two-particle op- 
erators with wave vectors P 7^ in the special case of 



iK^tP) : 



[d(P),H]; 
[dH-P),H] 



(158) 



ih-p(p) = ip(nn 

ih±D(P) = [D(P),H]; 

Following the equations of motion (158) we introduce 
four interconnected retarded Green's functions at T = 
[94, 95] 



Gii(P,t) = ((d(P,t);X+(P,0) / 
G 12 (P,t) = ((d\-P,t);X\P,0))); 
' 'p(P,t). 



G 13 (P,t) = 



-;Xt(P,0) 



(159) 



G 14 (At) = ((^;Xt (j P,0)\\; 



We also need their Fourier transforms Gij(P,u>), for 
which the equations of motion of the type similar to the 
equations of motion (158) were obtained. These Green's 
functions can be named as one operator Green's func- 
tions, because they contain only one two-particle oper- 
ator of the type d) , d, p, D. At the same time, in the 
right hand side of the corresponding equations of mo- 
tion there is a second generation of two-operator Green's 
function containing the different products of the two- 
particle operators mentioned above. For these operators 
was derived the second generation of the equations of 
motion containing in their right sides the Green's func- 
tion of the third generation. They are the three-operator 
Green's functions for which it is necessary to derive the 
third generation of equations of motion. However, we 
have to terminate here the evolution of the infinite chains 
of equations of motion for multi-operators Green's func- 
tion following the procedure proposed by Zubarev [95]. 
The truncation of the chains of the equations of motion 
and the decoupling of the one-operator Green's functions 
from the multi-operator Green's functions was achieved 
substituting the three operator Green's functions by the 
one-operator Green's functions multiplied by the aver- 
age value of remaining two operators. The average val- 
ues were calculated using the ground state wave function 
(153) and u-v transformations (155). The Zubarev proce- 
dure is equivalent to a perturbation theory with a small 
parameter of the type w 2 (l — v 2 ), which represent the 
product of a filling factor v — v 2 and the phase-space 
filling factor (1 — v 2 ) reflecting the Pauli exclusion prin- 
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ciple. The closed system of Dyson equations has the form 

4 

^G lj (P,u)E jk (P,u>) = C lk ; (160) 
k = 1,2,3,4 

There are 16 different components of the self-energy parts 
Ejfc(P, uj) forming a 4 x 4 matrix. Due to the structure 
of the self-energy parts the cumbersome dispersion equa- 
tion can be expressed in general form by the determinant 
equation 



det 12ij(P,u}) = 



(161) 



It splits into two independent equations. One of them 
concerns only the optical plasmon branch and has a sim- 
ple form 



E 33 (P,w) =0 



(162) 



It does not include the chemical potential jl and the 
quasiaverage constant rj. The second equation contains 
the self-energy parts En, E22, E44, E14, E41, E24 and 
E42, which include both parameters p, and fj. The sec- 
ond equation reads 



E 11 (P;w)E 2 2(P;w)E44( J P;^)- 
-E 4 i(P; w)E 22 (P; «)Su(P; w)- 
-E 4 2( J P;^)E 11 (P;w)E 2 4( J P;w) = 
The solution of the equation (162) is 

(MP)) 2 
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(D(P)D(-P)) 



N 2 

-U(Q-P))Sin 



Y,U{Q){U{-Q) - 

Q 
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The right hand side of this expression at small values of P 
has the dependence |P| 4 and tends to saturation at large 
values of P. The optical plasmon branch Hujop(P) has a 
quadratic dispersion law in the long wavelength limit and 
saturates in the rang e of short w avelengths. It depends 
on concentration as \fv 2 (\ — v 2 ) what coincides with the 



concentration dependencies of 3D plasma uj 2 



[96] and 2D plasma uj 2 (q) 



2 / ne 2 n 8 q 
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[88] , where n e and 



n s are the corresponding density of electrons. The sup- 
plementary factor (1 — v 2 ) in our case reflects the Pauli 
exclusion principle and the vanishing of the plasma os- 
cillations at v = v 2 = 1. The obtained dispersion law 
is shown in Fig. 2. Similar dispersion law was obtained 
for the case of 2D electron-hole liquid (EHL) in a strong 
perpendicular magnetic field [97], when the influence of 
the quantum vortices created by electron and hole sub- 
systems is compensated exactly. However, the saturation 



dependencies in these two cases are completely different. 
In the case of Bose-Einstein condensed magnetoexcitons 
it is determined by the ELLs, whereas in the case of EHL 
[97] it is determined by the Coulomb interaction in the 
frame of the LLLs. 

The acoustical plasmon branch has the dispersion law, 
which is completely different from the optical plasmon os- 
cillations. It has an absolute instability beginning with 
the small values of wave vector going on up to the con- 
siderable value pi w 2. In this range of wave vectors, the 
optical plasmons have energies which do not exceed the 
activation energy £7(0). It means that the optical plas- 
mons containing the opposite-phase oscillations of the 
electron and hole subsystems without displacement as 
a whole of their center of mass are allowed in the con- 
text of the attractive bath. On the other hand, the in- 
phase oscillations of the electron and hole subsystems in 
the composition of the acoustical plasmons are related 
to the displacements of their center of mass. Such dis- 
placements can take place only if their energy exceeds the 
activation energy U(P). As a result, the acoustical plas- 
mon branch has an imaginary part represented by the 
dashed line and is completely unstable in the region of 
wave vectors pi < 2. At greater values pi > 2 the energy 
spectrum is real and nonzero, approaching to the energy 
spectrum of the optical plasmons. 

In case of 2D magnetoexcitons in the BEC state with 
small wave vector kl < 0.5 described by the Hamiltonian 
(156), one should take into account that both continuous 
symmetries usual for the initial form (146) are lost. It 
happened due to the presence of the term J?(dt + rf-) in 
the frame of the Bogoliubov theory of the quasiaverages. 
Nevertheless the energy of the ground state as well as 
the self-energy parts Sjj(P, w) were calculated only in 

the simplest case of the condensate wave vector k = 0. 
These expressions can be relevant also for infinitesimal 



values of the modulus 



but with a well defined 



direction. In this case the symmetry of the ground state 
will be higher than that of the Hamiltonian (156), what 
coincides with the situation described by Georgi and 
Pais [33] . It is one possible explanation of the quasi-NG 
modes appearance in the case of exciton branches of the 
spectrum. Another possible mechanism of the gapped 
modes appearance is the existence of the local gauge 
symmetry, the breaking of which leads to the Higgs effect 
[7]. The interaction of the electrons with the attached 
vortices gives rise to a gapped energy spectrum of the 
collective elementary excitations as was established in 
Ref. [42, 54]. The number of the NG modes in the system 
with many broken continuous symmetries was deter- 
mined by the Nielsen and Chadha [17] theorem. It states 
that the number of the first-type NG modes Nj being 
accounted once and the number of the second type NG 
modes Nu being accounted twice equals or prevails the 
number of broken generators Nbg- It looks as follows 
Nj + 2Njj > Nbg- In our case the optical plasmon 
branch has the properties of the second-type NG modes. 
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We have Nj = 0; Nn = 1 and N B g = 2. It leads to the 
equality 2N H = N BG . The three branches of the energy 
spectrum are represented together in the Fig. 2. One 
of then is a second-type Nambu-Goldstone(NG) mode 
describing the optical plasmon-type excitations, the 
second branch is the first-type NG mode with absolute 
instability describing the acoustical-type excitations and 
the third branch is the quasi-NG mode describing the 
exciton-type collective elementary excitations of the 
system. 

We can repeat that results obtained in the magnetoex- 
citon system are similar to those obtained for the sys- 
tem of BEC of the quantum Hall excitons (QHExs) [81]. 
In these both models there is only one gapless Nambu- 
Goldstone mode between four branches of the energy 
spectrum. In our model it is related with the optical plas- 
mon branch, whereas in the case of QHExs this mode is 
represented by the superposition of the operators describ- 
ing the optical plasmon and exciton modes. In both mod- 
els the exciton branches of the spectrum are not gapless 
and differ from the NG modes. In our case the exciton en- 
ergy and quasienergy branches corresponding to normal 
and abnormal Green's functions have a gaps in the point 
k = 0, a roton-type segments in the region of interme- 
diary wave vectors kl ~ 1 and saturation-type behaviors 
at great values of kl. In the case of Ref[81] the exciton- 
type response function Xz(q,u) and the acoustical-type 
response function xf(q, w) have no poles in the region of 
small energies in the frame of the LLLs. It was concluded 
that the energies of these excitations may be situated at 
greater values. In our case the acoustical plasmon branch 
reveals an absolute instability in the range of small and 
intermediary values of k. It means that a real values of 
the pole does not exist in the range of small energies, 
which is similar with the results of Ref[81]. One can con- 
clude that the qualitative properties of the energy spec- 
tra in both models are similar in spite of the mentioned 
differences. It is an additional argument supporting the 
accuracy of our calculations, which satisfy to the Nielsen 
and Chadha theorem [17]. 

pi = 0; the second-type NG mode describing the opti- 
cal plasmons and the first-type NG mode with absolute 
instability (dotted line) describing the acoustical plas- 
mons. 

The result concerning the BEC at T = are estimates 
which describe the real situation at finite temperatures 
lower than the critical temperature of the Berezinskii- 
Kosterlitz-Thouless (BKT) topological phase transition 
[79, 80] related with the existence of the vortices and 
their clusters such as bound vortex-antivortex pairs. Just 
the unbinding of these pairs determines the critical tem- 
nnh 2 

perature Tbkt — -z — ; — , where n is the surface density 
2mkB 

of the Bose particles and m is their mass. On one side 
of the phase transition there is a quasi-ordered fluid and 
on the other is a disordered unbounded vortex plasma. 
Although the formation of an isolated vortex will not oc- 
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FIG. 2. Three branches of the collective elementary excita- 
tions: the exciton-type quasi-NG mode with a gap in the point 
pi = 0; the second-type NG mode describing the optical plas- 
mons and the first-type NG mode with absolute instability 
(dotted line) describing the acoustical plasmons. 

cur at low temperature due its extensive creation energy, 
there always can be production of a pair of vortices with 
equal and opposite charges since the perturbation pro- 
duced by such a pair falls off sufficiently rapidly at large 
distances so that their energy is finite [80]. Such topo- 
logical formations can be easily created by the thermal 
fluctuations. 

The presence of the vortex clusters makes the pre- 
viously infinite homogeneous 2D e-h system to become 
nonhomogeneous as a whole. However, the local homo- 
geneity with finite local surface areas can exist leading to 
the BEC with finite critical temperature T c = m kgig( n s) 
[98]. Instead of an off diagonal long-range order as in 
the case of 3D Bose gas in the 2D systems there is only 
a long rang correlations, which decays algebraically with 
distance. In such a way the quantum vortices promote 
the BEC and the formation of the superfluid component 
of the 2D Bose-gas at finite temperatures and at the same 
time the superfluid component is necessary for the for- 
mation of the quantum vortices. It is a self-organization- 
type situation. The BKT phase transition is a widely 
studied phenomenon [99-101]. 

Attempts to discover experimentally the spontaneous 
symmetry breaking in the exciton range of the spectrum 
and the efforts to evidence the spontaneous coherence in 
the 2D excitonic systems are will be considered in the 
next section on the basis of the references [102-146]. 

XII. SPONTANEOUS COHERENCE IN 2D 
EXCITONIC SYSTEMS 

As was mentioned by Snoke in Refs. [102, 103] re- 
cent experimental efforts of several groups have demon- 
strated the spontaneous coherence in polariton systems, 
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which can be viewed as a type of noncquilibrium BEC. 
The system of polaritons in the quantum wells embed- 
ded into the microcavity reveals the phenomenon of BEC 
and superfluity. The achievements in this field are pre- 
sented in Rcf. [104-110]. In these systems the polariton 
lifetime is longer than, but not much longer than the 
polariton-polariton scattering time, which leads to the 
thermalization. By contrast over past twenty years sev- 
eral groups of investigators represented by Snokc [111- 
121], Butov [122-130], Timofccv [131-138], Krivolapchuk 
[139-143], Fukuzawa [144-146] and their coworkers have 
pursued experiments in double quantum well (DQW) ex- 
citonic systems with very long lifetime. In these systems 
the indirect excitons (IXs) formed from spatially sepa- 
rated electrons and holes have dipole moments oriented 
perpendicularly to the layers. They are named dipole ex- 
citons and their interaction is not a short-range contact 
interaction but instead a long-range dipole-dipole repul- 
sion. We briefly recall the results obtained in Rcfs. [142, 
143]. 

When analyzing the possibility of BEC in a 2D sys- 
tem it should be noted that at T ^ condensation of 
a homogeneous 2D gas is impossible because of destruc- 
tion of the condensate by thermal fluctuations [147]. In 
a 2D system p(E) is constant and the integral N = 

CO 

J2 N k = J e(E -if/ksT-! would diverge at fj, ->• and 

k 

T ^ because of the zero denominator at the lower inte- 
grating limit and therefore BEC is impossible here. Phys- 
ically this fact means that the maximal occupation of free 
states (E > 0) is infinite. However, if in a 2D boson sys- 
tem, together with free excitons, there are present some 
discrete states (localized states whose existence is caused 
by the appearance of fluctuations of the heterointerface 
potential [148]) £o, £i etc, such that e < £i < E = 0, the 
situation changes essentially. In this case under increas- 
ing number of bosons in the system the value of chemical 
potential cannot be arbitrarily close to the E = value 
because of the N(eo) > requirement, so (— ^) m in = \so\ 
and, consequently the integral N has a finite value [149]: 

N 

n c (T) = ln(l - exp k B T ) ( i 65) 

Therefore, at the moment when n exceeds n c (T), local- 
ized states are occupied by a macroscopic number of par- 
ticles: n — n c (T) = n(eo) + n(si). This means that BEC 
into localized states occurs in a limited space region. In 
this sense BEC in a system of 2D bosons, having a dis- 
crete spectrum of energy together with a continuous one, 
resembles the experimentally discovered phenomenon of 
BEC in atoms of alkali metals in space-limited traps pro- 
duced by a magnetic field [150]. 

The obvious advantage of an IX in a DQW as a per- 
spective system to reveal BEC is the possibility of con- 
trolling effectively its radiative lifetime tr with the help 
of external effects. So, for example, an electric field Vdc 
applied to a DQW in the direction of the growth axis of 



the structure causes an essential decrease of the overlap 
of wavefunctions of an electron and a hole in an IX in 
the z direction and, as a result, tr increases significantly 
(by up to three orders of magnitude [145]). This allows 
a more effective cooling of the system to the bath tem- 
perature and, of equal importance, gives an opportunity 
to increase the concentration of the IX gas without in- 
creasing the pumping density. The latter circumstance 
plays an important role in the experiment since it allows 
us to decrease heating of the sample by phonons that are 
inevitably radiated at relaxation of photoexcited carri- 
ers and excitons. Just this heating of the sample under 
investigation is often the main cause of the impossibil- 
ity of reaching the critical temperature of the boson gas 
in experiments that use large optical pumping densities 
to create the critical density of bosons having very short 
lifetimes. 

Convincing evidence of the BEC effect would be an ap- 
propriate exciton distribution function over energies (mo- 
menta) obtained in an experiment. In general the exciton 
distribution function can be determined in experiment by 
the form of the phonon replica line in luminescence spec- 
tra, but in our study no phonon replica was observed. 
This was why in this paper the nonphonon luminescence 
line of space IXs was studied. However, since the inten- 
sity of exciton radiation is proportional to the occupation 
of radiative states by particles, it indirectly reflects the 
distribution function of excitons over the free and local- 
ized states, which both contribute to the formation of 
the inhomogeneously broadened IX line. Due to this fact 
one can hope that studies of IX luminescence will reveal 
the BEC effect predicted [149] for a system of 2D bosons 
which are distributed over the free and localized states. 

A giant (threefold) increase of luminescence intensity 
of a part of the spectral profile of the IX line in DQWs 
of GaAs/Al .33Gao.67As on changing the temperature of 
the sample and the value of the external electric field 
applied to DQWs was discovered. Besides that, the lu- 
minescence intensity of this part of the spectral profile 
of the IX line fluctuated with the characteristic time of 
tens of seconds. Such an unusual behaviour of the IX line 
was regarded as possible evidence for BEC in a system 
of 2D bosons placed in fluctuations of potential formed 
by heterointerfaces of the sample. 

Figure 1 of the Ref. [142] is reproduced in Fig. 3. It 
shows the luminescence spectra of the DQW, dependent 
on Vdc at T = 1.8K and the density of optical excitation 
P = hWcvnT 2 . Here, at Vdc = (figure (a)) the radiation 
spectrum was close to that of the flat-band case and it 
consisted of two lines, DXW and DXN, corresponding to 
luminescence of direct excitons (DXs) from the wide and 
the narrow wells respectively. At nonzero Vdc (figures 
(b)-(e)) an indirect regime was achieved (see the inset) 
when the IX line took the lowest energetic position in 
the PL spectra. When Vdc increased the IX line moved 
monotonically towards lower energies. 

It was noted that in some interval of Vdc a giant (up 
to threefold) increase (shot) of intensity of a part of 
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the IX spectral line profile (figure (d)) Ref. [142] oc- 
curred. A very important circumstance was that the 
intensity shot was absent from the whole investigated 
interval of spectral positions of the IX line at any tem- 
perature 4.2K < T < 30K and optical pumping densities 
of P < hWcrrT 2 . Thus the spectral profile of the IX line 
corresponding to the case of figure (d) was shown on a 
large scale in the figure 3 of Ref. [142] and it is repro- 
duced here in Fig. 4. It has some interesting peculiarities 
consisting of a narrow intense line C and of 'wings' W 
having significantly smaller intensity. Measurements of 
the temporal evolution of intensities of the C and W com- 
ponents have shown that the intensity of C, in contrast 
to that of W, fluctuated in time (changing threefold) on 
a characteristic scale of tens of seconds. 

The shape of the IX luminescence line in a DQW was 
inhomogeneous and was determined by exciton radiation 
from different space regions of the DQWs plane, which 
differed from one another by the thickness of the QW 
layers, by fluctuations of barrier composition and by the 
value of the local electric field of impurities. The emis- 
sion intensity of each spectral fragment of the IX line 
was proportional to the exciton occupation of the corre- 
sponding space region in the plane of the QW. However, 
a spectral region ( Vdc) where the situation changes es- 
sentially at T = 1.8K and P = hWcmT 2 appeared. In 
this region, shown in figure 3 (d) a significant (threefold) 
increase of luminescence intensity of a part of the spectral 
profile of the IX line and a consequent increase (1.5-fold) 
of integral intensity Ijx had been observed. Such a be- 
havior was anomalous in comparison to the monotonic 
decrease of Iix (and, most important, to the absence of 
the intensity shot) with increasing Vdc under other ex- 
perimental conditions (at T = 1.8K and P < lWcm~ 2 ; 
4.2K < T < 30K and any P < hWcmT 2 as well). This 
anomalous behavior of the IX line indicated that in the 
case of figure 1 (d) of Ref. [142] there are much more par- 
ticles participating in radiation (occupying states which 
can radiate) than in the cases of figures 1 (b), (c), (e) of 
Ref. [142]. 

To explain this the BEC model developed for 2D sys- 
tems [149] was proposed. It was shown that if in a system 
of bosons (excitons) there is a localized state £o below 
the bottom of the exciton band, the chemical potential 
of excitons is trapped by the localized level £o and, as a 
consequence, the number n c of particles appears to be fi- 
nite. Therefore just when the concentration n of excitons 
in the system under consideration exceeds n c , a macro- 
scopic number n — n c of particles comes to the lowest 
energetic state of the whole boson system (i.e. free zone 
and the localized state) and that leads to the appearance 
of a condensate. 

The phenomenon of a giant intensity shot was revealed 
in an experiment when the IX line shifted to lower ener- 
gies (that means changing Vdc), this being a consequence 
of an increase of the radiative lifetime of excitons due 
both to an increase of the exciton concentration (at con- 
stant optical pumping) and to their effective thermaliza- 
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FIG. 3. PL spectra taken at T = 1.8K, P = 5Wcm' 2 and 
V dc = V(a), -0.5 V(b), -2 V(c). The inset shows the indirect 
regime of the DQW following the Ref. [142]. 
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FIG. 4. IX line spectral profile of figure (d). Two vertical dot- 
ted lines separate different spectral parts of the IX line profile 
which shown (C) and did not show (W) temporal evolution 
of the PL intensity. The data were obtained in Ref. [142] . 

tion to the bath temperature. These two circumstances 
according to the authors opinion cause the BEC that 
leaded to a huge population of a localized state. 
Since, as noted above, all localized excitons take part 
in radiative recombination, this resulted in a significant 
change of the shape of the IX luminescence line from the 
case of figure (c) to that of figure (d) . Thus the totality 
of the experimental data describing the evolution of the 
IX luminescence line (namely, a giant rise of intensity of 
a part of the spectral profile of the IX line followed by 
long-time oscillations) in the authors opinion provided 
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evidence that in a system of IXs of high density a Bosc- 
Einstcin condensate at localized states (traps) in DQWs 
appeared. Thermal equilibrium of this type of excitons in 
a trap has been demonstrated experimentally [102, 103]. 

As was appreciated in the Ref. [102, 103] up to now 
there has not been an universally accepted demonstration 
of BEC in this type of systems and is necessary a better 
understanding of the many-body effects of the interact- 
ing dipole IXs. But the accumulated knowledge permits 
to formulate some conclusions. One of them states that 
the confinement of the excitons in a trap analogous to 
the optical traps for the cold atoms is a great advantage 
instead of creating excitons with a laser allowing them 
to expand freely out of the excitation region. Recent 
work [117] showed that the IXs in DQWs reach equilib- 
rium both energetically and spatially in a stress induced 
trap. One variant of the BEC of excitons equivalent to 
the BEC in a trap was proposed by Jan and Lee [149] 
and was used in Ref. [142, 143]. Another conclusion 
formulated in Ref. [102. 103] concerns the role of the 
temperature. If the temperature is low compared to the 
energy fluctuation due to the disorder then the excitons 
will become trapped in low energy minima of the disor- 
der potential and will not act as a free gas. Such energy 
minima can localize only one or a small number of IXs 
because they are repealing each other due to the dipolc- 
dipole interaction. In difference on it in the trap there 
are localized levels able to accommodate a macroscopic 
number of IXs. A strong energy shift due to interactions 
may cancel out the trapping potential and may flatten it 
[121]. 

Another important conclusion of Rcfs. [102, 103] is 
based on the results by Laikhtman and Rapaport [151] 
who underlined that the dipole IXs in a coupled quan- 
tum wells (CQWs) no longer act as a gas but rather as 
a correlated liquid. This does not mean that BEC is im- 
possible at high density. But the canonical telltale for 
condensation in a weakly interacting Bose gas, namely 
a peak of occupation number at k = may not be eas- 
ily seen in these systems. It would be better to look for 
hydrodynamic effects of condensation of excitons such as 
quantum vortices or superfluidity similar to the liquid 
helium. 

As one can see the most eminent achievements dis- 
cussed in the review article concern the famous FQHEs 
discovered in the frame of the one-component 2DEG. 
They suggest to search similar phenomena in the frame of 
the two-component, 2D e-h systems, when the CPs will 
be formed by electrons and holes with attached point vor- 
tices in different combinations taking into account the in- 
teractions inside the CPs from one side and the Coulomb 
interaction between the electrons and holes forming the 
usual magnetoexcitons from the other side. In spite of 
the fact that the supplementary gauge magnetic fields 
created by electrons and holes with opposite charges will 
be compensated in the mean field approximation, never- 
theless the new electron quantum states will appear as 
elementary excitations and quantum fluctuations in these 



new conditions. 



XIII. CONCLUSIONS: 

The purpose of the present review is to discuss the 
phenomena related to the spontaneous breaking due to 
the quantum fluctuations of the continuous symmetries 
existing in the frame of the two-dimensional e-h systems 
in a strong perpendicular magnetic field with electrons 
and holes lying on the lowest Landau levels. The spon- 
taneous symmetry breaking leads to the formation of the 
new ground states and phase transitions and determines 
the energy spectra of the collective elementary excita- 
tions appearing over the new ground states. 

The main attention is given to the electron-hole sys- 
tems forming the coplanar magnetoexcitons in the Bose- 
Einstcin condensation ground state with wave vector 
k = under the influence of the excited Landau lev- 
els when the exciton-type excitations coexist with the 
plasmon-type oscillations. At the same time the prop- 
erties of the 2DEG under the conditions of the FQHE 
as well as of the similar 2DHG spatially separated on 
the layers of the DQW are taken into account, so as 
to foresee their possible influence on the states of the 
coplanar magnetoexcitons when the distance between the 
DQW layers diminishes. Side by side with the 2DEG 
and 2DHG a bilayer electron systems in the conditions 
of the FQHE with one half filling factor of LLLs in each 
layer and with the total filling factor of two layers equal 
to unity are taken into account because the coherent 
superposition of the electron states in two layers hap- 
pens to be equivalent with the formation of the QHExs 
in the coherent macroscopical state, which can be com- 
pared with the BEC of the coplanar magnetoexcitons. 
The breaking of the global gauge symmetry as well as of 
the continuous rotational symmetries leads to the forma- 
tion of the gapless Nambu-Goldstone modes of the col- 
lective excitations above the selected ground state, corre- 
sponding to the macroscopical wave function with a fixed 
phase, whereas the breaking of the local gauge symmetry 
gives rise to the Higgs phenomenon characterized by the 
gapped branches of the energy spectrum of the collective 
elementary excitations. The existence of the gapless and 
gapped branches of the energy spectrum is equivalent to 
the appearance of the massless and massive particles cor- 
respondingly in the relativistic physics. 

Application of the Nielsen-Chadha theorem establish- 
ing the relation between the number of the NG modes 
and the number of the broken symmetry operators as well 
as the elucidation of the conditions when the quasi-NG 
modes appear was effectuated on the concrete example 
of the spinor atoms in the state of BEC in an optical trap 
with the aim to better understand the results concerning 
the coplanar magnetoexcitons. The Higgs phenomenon 
gives rise to the formation of the composite particles in 
the frame of the 2DEG in conditions of the FQHE, so 
that the electron with an odd or even number of the at- 
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tached point vortices behaves as an composite boson or 
fcrmion correspondingly. Their description in the frame 
of the Ginzburg-Landau theory is demonstrated. 

Side by side with the 2D coplanar magnetoexcitons 
the conditions under which the spontaneous coherence 
may appear in the system of indirect excitations in the 
double quantum well structures with spatially separated 



electrons and holes were discussed. The experimental at- 
tempts to achieve the BEC of IXs in the traps arising due 
to the interface width fluctuations or due to the applied 
stress were reviewed, the concluding remarks and the 
recommendations are mentioned. The formation of the 
high density 2D magnetoexcitons and magnetoexciton- 
polaritons with point quantum vortices attached is sug- 
gested. 
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